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Abstract 

Comparative  analysis  is  the  problem  of  predicting  how  a  system  will  react 
to  perturbations  in  its  parameters,  and  why.  For  example,  comparative  analysis 
could  be  asked  to  explain  why  the  period  of  an  oscillating  spring/block  system 
would  increase  if  the  mass  of  the  block  were  larger.  This  paper  formalizes  the 
problem  of  comparative  analysis  and  presents  a  technique,  differential  qualita¬ 
tive  (DQ)  analysis,  which  solves  the  task,  providing  explanations  suitable  for 
use  by  design  systems,  automated  diagnosis,  intelligent  tutoring  systems,  and 
explanation  based  generalization. 

DQ  analysis  uses  inference  rules  to  deduce  qualitative  information  about  the 
relative  change  of  system  parameters.  Multiple  perspectives  are  used  to  repre¬ 
sent  relative  change  values  over  intervals  of  time.  Differential  analysis  has  been 
implemented,  tested  on  a  dozen  examples,  and  proven  sound.  Unfortunately, 
the  technique  is  incomplete;  it  always  terminates,  but  does  not  always  return 
an  answer. 
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Comparative  analysis  is  the  problem  of  predicting  how  a  system  will  react 
to  perturbations  in  its  parameters,  and  why.  For  example,  comparative  analysis 
could  be  asked  to  explain  why  the  period  of  an  oscillating  spring/block  system 
would  increase  if  the  mass  of  the  block  were  larger.  This  paper  formalizes  the 
problem  of  comparative  analysis  and  presents  a  technique,  differential  qualita- 
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tive  (DQ)  analysis,  which  solves  the  task,  providing  explanations  suitable  for 
use  by  design  systems,  automated  diagnosis,  intelligent  tutoring  systems,  and 
explanation  based  generalization. 

DQ  analysis  uses  inference  rules  to  deduce  qualitative  information  about  the 
relative  change  of  system  parameters.  Multiple  perspectives  are  used  to  repre¬ 
sent  relative  change  values  over  intervals  of  time.  Differential  analysis  has  been 
implemented,  tested  on  a  dozen  examples,  and  proven  sound.  Unfortunately, 
the  technique  is  incomplete;  it  always  terminates,  but  does  not  always  return 
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1  Introduction 


The  problem  of  symbolic  analysis  of  real-world  systems  is  central  to  many  prob¬ 
lems  in  artificial  intelligence.  In  order  to  cope  with  a  changing  world  one  must 
be  able  to  understand  its  behavior.  Recently,  considerable  emphasis  has  been 
put  on  a  specific  kind  of  analysis:  qualitative  simulat  ion[2.7.2fi.  1 5.21].  Qual¬ 
itative  simulation  seeks  to  produce  a  description  of  t In-  behavior  of  a  system 
over  time,  often  in  the  form  of  a  tree  of  histories  of  the  system's  qualitatively 
interesting  changes  over  time  [23], 

This  paper  discusses  the  problem  of  comparative  analysis,  in  many  ways 
the  complement  of  qualitative  simulation,  and  describes  an  implemented,  sound 
solution  technique  called  differential  qualitative  (DQ)  analysis  Whereas  qual¬ 
itative  simulation  takes  a  structural  description  of  a  system  and  predicts  its 
behavior,  comparative  analysis  takes  as  input  this  behavior  and  a  perturbation 
and  outputs  a  description  of  how  and  why  the  behavior  would  change  as  a  result 
of  the  pertubation. 

For  example,  given  the  structural  description  of  a  horizontal,  frictionless 
spring/block  system  (e  g.,  Hooke's  law),  a  qualitative  simulator  would  say  that 
the  block  would  first  move  one  direction,  then  stop,  then  reverse,  etc.  A  de¬ 
scription  of  oscillation  would  result.  Comparative  analysis,  on  the  other  hand, 
takes  this  description  of  oscillation  and  evaluates  the  effects  of  perturbations. 
For  example,  it  would  deduce  that  the  period  would  lengthen  if  the  mass  of 
the  block  were  increased,  and  explain  why.  Just  as  qualitative  simulation  works 
without  explicit  equations  for  the  value  of  each  parameter  as  a  function  of  time, 
comparative  analysis  does  not  need  a  formula  for  the  period  of  oscillation. 

The  importance  of  the  qualitative  approach  to  comparative  analysis  is  the 
resulting  explanation  of  why  the  behavior  changes  If  it  weren't  for  the  explana¬ 
tion,  one  might  simply  solve  a  differential  equation  model  using  using  symbolic 
or  numeric  techniques.  Many  artificial  intelligence  problems,  for  example  de¬ 
sign,  diagnosis,  and  intelligent  tutoring  systems,  have  comparative  analysis  as 
an  important  component;  the  explanation  is  used  in  many  different  ways. 

•  One  way  method  of  automated  design  is  the  principled  modification  of 
previous  designs  [25].  For  example,  suppose  a  library  design  for  a  VLSI 
pullup  circuit  has  too  long  a  rise  time.  If  the  problem  solver  considers 
increasing  the  width  of  some  wire  to  decrease  the  rise  time,  it  would  like 
to  know  the  ramifications  of  this  modification  relative  to  the  initial  be¬ 
havior.  Will  the  delay  decrease?  What  happens  to  power  dissipation? 
Comparative  analysis  answers  these  questions  in  qualitative  terms,  as  is 
appropriate  for  initial  design  evaluation.  Hy  analyzing  an  explanation 
for  why  the  changes  happen,  the  problem  solver  could  focus  on  further 
changes  to  counteract  undesired  effects 

•  M  any  of  the  programs  which  perform  diagnosis  from  first  principles  use 
similar  generate  and  test  paradigms  [!•].  Comparative  analysis  ran  simplify 


diagnosis  of  continuous  systems  (such  as  analog  electronics)  in  two  ways. 
Comparative  analysis  provides  a  direct  test  for  certain  hypothesized  faults; 
if  one  suspects  a  resistor  of  a  low  value,  comparative  analysis  can  predict 
the  resulting  behavior.  If  this  prediction  does  not  match  the  observed 
behavior,  the  generator  might  use  the  explanation  to  suggest  or  rule  out 
additional  candidate  faults. 

In  addition,  the  specific  type  of  comparative  analysis  discussed  in  this 
paper,  DQ  analysis,  can  be  used  backwards  to  generate  candidate  faults. 
If  an  output  voltage  measures  too  low.  reversing  the  inference  rules  of 
section  3  might  lead  to  the  hypothesis  that  some  capacitor  has  too  high  a 
value. 

•  A  key  subproblem  of  intelligent  tutoring  systems  (ITS)  is  the  automatic 
explanation  of  the  behavior  of  complex  systems.  Most  AI  work  in  this  di¬ 
rection  has  focused  on  the  role  of  qualitative  simulation  when  explaining 
the  mechanism  through  which  devices  achieve  functionality  [19,8].  Qual¬ 
itative  simulation  is  a  critical  component  of  explanation  generation,  but 
understanding  how  systems  respond  to  changes  is  also  important.  One 
doesn't  really  understand  the  workings  of  a  refrigerator,  if  one  can’t  ex¬ 
plain  the  effect  of  a  stronger  compressor  on  efficiency  and  minimum  tem¬ 
perature. 

The  rest  of  this  paper  shows  how  L)Q  analysis  can  solve  comparative  analysis 
problems  and  produce  clear  explanations  as  well.  The  trick  to  DQ  analysis 
is  the  use  of  multiple  ‘perspectives'  to  define  relative  change.  Inference  rules 
manipulate  these  relat  ive  change  values  to  generate  causal  arguments  that  solve 
comparative  analysis  questions 

The  rest  of  the  introduction  explains  how  differential  qualitative  analysis 
solves  comparative  analysis  problems  Section  1.1  presents  more  detail  about 
the  spring/block  example1  to  illustrate  the  important  notion  of  ‘perspective’. 
Section  1.2  introduces  a  heat  exchanger  example  to  emphasize  the  importance 
of  considering  multiple  behavioral  topologies.  Together  these  two  sections  show 
the  range  of  questions  that  the  differential  approach  to  comparative  analysis 
can  answer.  Section  1.3  suggests  a  different  approach  to  comparative  analysis: 
a  novel  technique  called  exaggeration  Finally,  section  1.4  gives  an  overview  of 
the  remainder  of  the  paper 

1.1  Perspectives 

Perspectives  are  the  most  important  concept  in  DQ  analysis:  they  are  best 
introduced  with  an  example  Consider  an  ideal  spring  attached  to  a  block  on  a 
frictionless  table  (figure  1) 


1  Alt  the  examples  in  this  paper,  an. I  a  d.r/en  mere,  have  heeri  implemented  and  tested. 


I  fl  II  i  I  FTTITI  M  Mfl"l  I1  II  M  >  I  JI.I.i.i.ii.  i 


fUiUWUWtUfUi^yrJF)WWlW!Wipi|piW^fMI»B,PmnJfj 


Figure  1:  Ideal  Spring  Attached  to  Block  on  Frictionless  Table 

The  system  can  be  defined  in  Kuipers’  QSIM  [15]  notation  in  terms  of  six 
parameters,  each  a  function  of  time:  spring  constant  —  A’,  mass  M,  position 
A',  velocity  V,  acceleration  .4,  and  force  F  related  by  Newton's  second  law 
{F  =  M A)  and  Hooke’s  law  (F  =  —KX).  Mass  and  spring  constant  are 
independent  parameters  that  remain  at  constant  values  over  time.  The  initial 
conditions  are  specified  as  follows:  A/(0)  >  0,  —  A'(0)  <  0,  V'(0)  =  0,  and 
Af(0)  =  x0  <  0. 

This  description  may  now  be  simulated,  but  because  of  ambiguities  inher¬ 
ent  in  qualitative  values  [14],  QSIM  produces  several  possible  behaviors  for  this 
system,  including  ones  corresponding  to  increasing,  decreasing,  and  stable  oscil¬ 
lation.  Although  comparative  analysis  could  be  done  on  all  of  these  behaviors, 
for  this  example,  I  assume  the  interpretation  of  stable  oscillation  (figure  2). 

Now  we  are  ready  to  pose  a  comparative  analysis  problem. 

Example  1  What  happens  to  the  period  of  oscillation  if  the  mass  of  the  block 
is  increased? 

The  answer  is  that  the  length  of  the  period  increases: 

Since  force  is  inversely  proportional  to  position,  the  force  on  the 
block  will  remain  the  same  when  the  mass  is  increased.  But  if  the 
block  is  heavier,  then  it  won’t  accelerate  as  fast.  And  if  it  doesn't 
accelerate  aa  fast,  then  it  will  always  be  going  slower  and  so  will 
take  longer  to  complete  a  full  period  (assuming  it  travels  the  same 
distance). 

What  kind  of  information  is  needed  to  produce  this  explanation7  Take  the 
first  step:  “The  force  on  the  block  will  remain  the  same.”  Figure  3  shows  a  real¬ 
valued  plot  of  force  versus  time.  The  graph  of  force  in  the  perturbed  system  is 
drawn  with  a  dotted  line;  in^the  text  1  will  distinguish  the  two  parameters  by 
calling  the  perturbed  force  F. 

Clearly,  F  F  as  a  function  of  time.  The  corresponding  values  of  F  and 
F  are  different  for  almost  every  possible  time.  The  real  meaning  of  “The  force 
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Figure  2:  QSIM  Behavior  for  Stable  Spring  Oscillation 


on  the  block  will  remain  the  same”  is  that  F  and  F  are  the  same  for  all  values 
of  X.  Although  this  reparametrization  was  not  mentioned  explicitly  in  the 
explanation,  it  is  essential  to  the  soundness  of  the  argument. 

In  order  to  allow  programs  to  generate  and  evaluate  explanations  like  the 
one  for  the  spring  and  block,  it  is  necessary  to  take  this  implicit  concept  and 
make  it  explicit.  I  do  this  with  the  use  of  ‘perspectives’.  Thus  the  first  line  of 
the  argument  could  be  rewritten  “If  the  mass  is  increased,  force  does  not  change 
from  the  perspective  of  position.”  Making  perspectives  explicit  is  the  crucial 
step  in  performing  DQ  analysis  to  solve  a  comparative  analysis  problem.  Once 
the  notion  of  perspective  is  explicit,  one  can  address  questions  like  “Which  per¬ 
spective  best  suits  a  problem”  and  “What  inferences  are  sound?”  The  answers 
are  not  as  obvious  as  they  might  appear. 

For  example,  consider  the  ‘obvious’  inference  “Since  it  is  going  slower  it  will 
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Figure  3:  Actual  Plot  of  Force  versus  Time 


take  longer  to  go  the  same  distance."  But  what  does  it  mean  for  the  block  to  be 
going  slower?  From  what  perspective  is  velocity  lower?  If  velocity  were  lower 
from  the  perspective  of  time,  then  the  conclusion  would  indeed  be  obvious.  But 
just  as  with  the  parameter  force  (figure  3),  there  are  times  when  the  perturbed 
velocity  is  not  lower  than  it  was  in  the  original  system.  Once  again,  position 
is  the  correct  perspective.  In  fact,  as  shown  in  section  3.  the  explanation  is 
correct,  but  it  would  not  necessarily  be  so  if  the  perspective  was  some  other 
parameter. 

Reasoning  about  perspectives  explicitly,  and  using  sound  rules  of  DQ  analy¬ 
sis  (section  3),  the  CA  program  has  correctly  generated  the  correct  solution  and 
an  explanation  like  that  shown  above.  Here  is  another  example  which  it  solves 
by  using  perspectives  in  a  different  way. 


Example  2  H/ial  happens  to  the  maximum  velocity  if  the  initial  displacement 
is  increased? 

C A  generates  the  justification  which  can  be  turned  into  English  as  follows: 

Since  K  and  M  haven’t  changed,  the  force  on  the  block  is  the 
same  for  any  position  that  the  block  used  to  pass  through.  So  the 
acceleration  is  the  same  for  any  position.  But  since  the  initial  dis¬ 
placement  has  been  increased,  the  block  will  already  be  moving  when 
it  reaches  the  old  initial  position,  where  previously  the  block  was 
stopped.  Since  the  accelerations  are  the  same  from  here  on,  and  the 
block  is  already  moving  faster,  it  will  keep  on  moving  faster  and  will 
have  a  higher  maximum  velocity. 

The  rules  which  compose  this  reasoning  are  explained  in  section  3. 

1.2  Changes  in  Behavioral  Topology 

The  previous  section  showed  how  the  explicit  use  of  perspectives  could  determine 
the  relative  change  of  parameter  values  and  time  durations  given  an  initial 
perturbation.  However,  sometimes  the  perturbation  results  in  change  of  a  more 
fundamental  nature.  Consider  the  heat  exchanger  shown  in  4.  Hot  oil  flows 
through  the  pipe  losing  heat  to  the  cold  water  bath  as  it  goes.  Figure  5  shows 
a  possible  QSIM  behavior  that  corresponds  to  the  rase  when  the  hot  oil  reaches 
thermal  equilibrium  just  as  it  exits  from  the  pipe.  (Remember  that  since  this 
is  a  qualitative  plot,  the  apparent  slope  does  not  imply  that  these  functions  are 
linear.)  Let’s  pose  a  comparative  analysis  problem. 


Figure  4:  Hot  Oil  Flows  Through  Heat  Exchanger 
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Figure  5:  Initial  Behavior  of  Heat  Exchanger 


Example  3  What  happens  to  the  behavtor  of  the  heat  exchanger  if  the  thermal 
conductivity  is  increased ? 

The  answer  is  that  the  oil  will  more  more  quickly  than  before.  And  since 
the  oil  is  flowing  through  the  exchanger  at  the  same  rate,  it  must  reach  ther¬ 
mal  equilibrium  before  leaving  the  pipe  (figure  6).  Thus,  unlike  the  previous 
examples  where  the  perturbation  resulted  only  in  continuous  changes  in  various 
parameters,  the  perturbation  of  example  3  caused  a  discontinuous  change:  the 
previously  cotemporaneous  ‘events’  of  thermal  equilibrium  and  disgorgement 
from  the  pipe  now  happen  at  different  times. 

I  call  the  switch  from  figure  5  to  6  a  change  in  behavioral  topology.  Example 
3  is  a  simple  case  of  topological  change:  the  initial  behavior  was  inconsistent  and 
a  single  new  behavior  was  indicated.  However,  the  situation  isn't  always  so  easy. 
Section  4  describes  how  perturbations  can  lead  to  multiple  consistent  behaviors 
and  presents  heuristics  for  determining  the  most  likely  resulting  behavior 

1.3  Exaggeration 


While  most  of  this  paper  deals  with  the  DQ  solution  technique  to  comparative 
analysis  problems,  it  is  worth  noting  that  other  qualitative  techniques  can  solve 
similar  problems  One  such  technique,  called  exaggeration  [20],  produces  expla¬ 
nations  that  are  completely  different  from  those  of  the  differential  technique 
Consider  the  question  of  example  1:  “What  happens  to  the  period  of  oscillation 
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Figure  6:  Perturbed  Behavior  of  Heat  Exchanger 

if  the  mass  of  the  block  is  increased?"  Compare  the  exaggerated  explanation 
with  the  one  generated  by  DQ  analysis. 

If  the  mass  were  infinite,  then  the  block  would  hardly  move  at 
all.  So  the  period  would  be  infinite.  Thus  if  the  mass  was  increased 
a  bit,  the  period  would  increase  as  well. 

Exaggeration  is  a  kind  of  asymptotic  analysis — the  perturbation  is  taken  to 
the  limit  to  make  the  effect  more  easily  visible.  Exaggeration  is  common  in 
intuitive  descriptions  of  physical  behavior  and  appears  quite  powerful.  As  the 
example  shows  it  often  results  in  a  concise  explanation. 

But  exaggeration  is  subtle.  It  works  only  when  the  system  responds  mono- 
tonically  to  perturbations.  Furtnermore,  it  requires  non-standard  analysis  to 
reason  about  infinity.  It  s  quite  easy  to  concoct  a  plausible  exaggerated  argu¬ 
ment  which  is  faulty,  and  a  careful  formalization  of  the  technique  is  beyond  the 
scope  of  this  paper  See  [20,22]  for  details. 

1.4  Overview 

The  next  section  is  foundational — it  shows  how  perspectives  are  essential  to  a 
meaningful  definition  of  relative  change  Section  3  explains  how  the  differential 
approach  to  comparative  analysis  can  be  implemented  by  a  number  of  inference 
rules  The  rules  are  proved  sound,  and  their  adequacy  is  discussed.  Section 


1  slews  how  to  j  <  r«'>  1  i> '  t  tin1  effect  of  perl  urhat  ions  that  change  the  behavioral 
topology  of  a  system  In  section  a  i  discuss  the  relationship  between  com- 
parati\e  analysis  and  previous  work  in  mathematics,  engineering,  and  artificial 
intelligence  the  paper  concludes  with  a  discussion  of  areas  for  future  research 
til  l  is  followed  by  a  brief  appendix  that  contains  details  too  technical  for  the 
main  b.  m|\  of  t  lie  text 


2  Preliminaries 


As  my  formalism  is  based  on  that  used  by  Kuipers  for  QS1.M  [15],  I  start  out 
by  summarizing  his  definitions. 

Definition  1  .4  PARAMETER  is  a  reasonable  function  of  time. 

See  [15]  for  the  actual  definition  of  reasonable  function;  the  intuition  is  that 
of  continuity,  continuous  differentiability,  and  a  finite  number  of  critical  points 
(places  where  its  derivative  is  zero)  Parameters  are  denoted  by  capital  letters. 
Thus  the  velocity  of  a  projectile  might  be  described  by  the  parameter,  V,  which 
is  a  function  that  maps  time  to  velocity. 

Definition  2  Each  parameter  has  an  associated  set  o/LAN'DMARK  VALUES  which 
is  a  subset  of  the  range  of  the  parameter.  The  landmark  values  always  include 
(but  aren  't  restricted  to)  zero,  the  values  of  the  parameter  at  the  beginning  and 
ending  times ,  and  the  values  of  the  parameter  at  each  of  its  critical  points.  A 
tune,  t,  is  a  DISTINGUISHED  TIME  POINT  of  a  parameter  P  if  it  ts  a  boundary 
element  of  the  set  of  tunes  that  P(t)  =  p ,  for  some  landmark  value 

Landmark  values  are  those  values  considered  to  be  interesting  to  the  human 
observer,  and  the  times  when  these  values  are  reached  are  of  interest  too.  When 
a  parameter  becomes  constant  for  an  interval  of  time,  then  it  will  take  on  a 
landmark  value  for  infinite  number  of  time  points.  This  is  why  the  definition 
only  considers  the  boundary  times  distinguished. 

Definition  3  A  SYSTEM  is  a  set  of  parameters  that  are  related  with  a  STRUC¬ 
TURAL  DESCRIPTION  that  consists  of  a  finite  set  of  qualitative  differential  equa¬ 
tions  defined  using  the.  following:  turn  differentiation,  addition,  multiplication, 
and  relation  by  monotonia  functions. 

Kuipers’  program.  QSIM.  takes  a  system  and  a  set  of  initial  values  for  each 
of  the  parameters  and  produces  a  set  of  possible  behaviors  for  the  system;  the 
definitions  below  describe  this  behavioral  output: 

2.1  Qualitative  Behavior 

Definition  4  Let  pa  <  <  pe  be  the  landmark  values  of  a  parameter  P.  For 

any  time  t  define  the  value  of  /’  at  t  as: 

QYA  UP.,)  =  l  ,  ,f  Pit)  =  landmark  Vl 

\  (P)  P}+ i )  ij  i  (i)  e  (p}.p,  + 1 ) 

Define  the  direction  of  P  at  I  as: 
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inc 

if  7i  P(t)  >  0 

QDIR(  P,t)  —  std 

>f&P(t)=0 

dec 

if  It  P(l)  <  0 

Define,  QS(P.t),  the  state  of  P  at  t,  as  the  pair:  <QVAL(P,t).  QDIR(P,t)> 

The  qualitative  state  over  the  interval  between  two  adjacent  distinguished 
time  points  is  defined  similarly. 

Defi  nition  5  For  any  parameter  P ,  the  behavior  of  P  is  a  sequence  of  states 
of  P : 

qs(P,  t0),  QS(P,  t0,  <i),  QS(P,ti) . QS(P,tn-i'tn),QS(P.tn) 

alternating  between  states  at  distinguished  time-points,  and  states  on  intervals 
between  distinguished  time-points. 

Recall  that  a  system  contains  a  set  of  parameters  each  with  its  own  land 
marks  and  distinguished  time  points. 

Definition  6  The  DISTINGUISHED  TIME-POINTS  of  a  system  are  the  union  of 
the  distinguished  time-points  of  the  parameters.  Thus  the  state  of  a  system 
changes  whenever  the  state  of  any  parameter  changes.  The  BEHAVIOR  of  a  sys¬ 
tem  is  thus  a  sequence  of  system-states  alternating  between  distinguished  time- 
points  and  intervals. 

To  perform  comparative  analysis  it  is  necessary  to  abstract  away  from  spe¬ 
cific  times,  since  two  different  systems  may  have  analogous  behaviors,  but  change 
states  at  different  times.  This  is  where  my  formal  treatment  diverges  from  that 
of  Kuipers. 

Definition  7  A  parameter  is  said  to  reach  a  TRANSITION  when  its  QVAL  changes 
to  or  from  a  landmark  value.  .4  system  ts  said  to  reach  a  TRANSITION  when  any 
parameter  transitions.  Transitions  only  occur  at  distinguished  lime-points,  and 
every  distinguished  time  point  marks  a  transition.  It  will  prove  useful  to  be 
able  to  refer  to  these  transitions  independent  of  the  time  at  which  they  occur, 
thus  the  sequence  of  transitions  for  a  behavior  will  be  denoted  by  the  set  {-,} 
Every  behavior  also  has  a  TIME  FUNCTION.  T .  which  takes  transitions  tv  the 
distinguished  lime-points  when  they  occur. 

The  intuition  is  that  each  7  marks  an  event  which  changes  the  state  of  the 
system.  When  comparing  two  behaviors.  I  match  them  up  event  by  event  and 
use  the  time  functions  to  tell  whether  one  system  is  changing  faster  or  slower 
than  the  other. 


■  -  ‘ 


2.2  Comparing  Two  Beha\  iors 

To  compare  two  behaviors,  they  must  be  distinguishable:  I  use  the  hat  accent 
to  denote  the  second  behavior.  Thus  T  denotes  the  time  function  of  the  second 
system,  and  F{T{~n  ))  denotes  the  second  system’s  value  of  F  at  the  time  of  the 
first  transition  To  simplify  the  problem  of  of  comparative  analysis,  I  -tart  by 
only  comparing  systems  with  identical  structural  descriptions  whose  behaviors 
are  topologically  equal,  as  defined  below 

Definition  8  The  behaviors  of  tuo  system*.  S  and  S.  are  TOPOLOGICALLY 

EQUAL  if  they  have  the  same  sequence  of  transitions.  70 . *■*..  and  foratl  1 

such  that  0  <  i  <  k . 

qs(5.T(7,)1  -  qs(S, T (*,)) 

and  forall  i  such  that  0  <  i  <  Ic, 

qs(.>'.  T'7i).T(7,  +  1))  -  qs(5.  f  \  ;  i-1+1)) 

The  assumption  of  topological  equality  rules  out  possibilities  like  the  block 
failing  to  make  a  complete  oscillation  if  its  mass  was  increased  too  much,  but 
it  does  allow  a  certain  pliability.  If  two  behaviors  are  topologically  equal,  their 
respective  sets  of  landmarks  share  the  same  ordinal  relationships,  but  the  un¬ 
derlying  real  values  for  the  landmarks  can  be  different. 

Section  4  explains  how  this  assumption  can  be  relaxed,  but  even  with  it,  the 
problem  is  nontrivial  Consider  two  oscillating  spring  block  systems.  Even  if  the 
blocks  have  different  mass  and  the  spring  constants  differ,  the  two  systems  have 
topologically  equal  behavior.  Yet  the  relative  values  of  parameters  such  as  period 
of  oscillation  may  be  different  These  are  the  first  changes  that  comparative 
analysis  must  determine. 

Before  I  can  explain  the  techniques  for  performing  comparative  analysis,  I 
need  to  present  a  notation  for  describing  t h« •  desired  output  Its  easy  to  compare 
the  values  of  parameters  at  transition  points 

Definition  9  Given  a  parameter.  F  an  !  a  transition  define  the  RELATIVE 
CHANGE  (RC)  of  l  at  7,  as  follows 

fh,  tf\fm->,))\  >  |/'(T(7.))i 

n\,  >f\F(Th,))\  =  \Frr(y,))\ 

Fll,  tf\F(T(-n))\  <  \F(T( -„))( 

For  example,  if  the  t  wo  spring  bh .ek  s\ stems  w,  n  I ,ot h  st ai  led  with  negative 
displacement  and  zero  velocity  (i.e.  A  <  C  and  l  =  0).  th.-ir  first  transition 
would  occur  when  A  readied  zero  I  his  notation  allows  one  to  express  that  the 
second  block  is  moving  slower  at  the  p  in!  >>f  1  ran»it ion  IT1;  It  is  important 
to  distinguish  the  relative  change  notation  troiii  statements  about  values  and 


derivatives.  E.'ven  though  QVALfl',  T(-)i))  is  positive,  and  QDIR(V',  T(-)i )) 
is  std. 

The  curious  reader  may  wonder  at  the  use  of  absolute  values  in  this  defini¬ 
tion.  Relative  change  could  also  be  defined  by  comparing  signed  values.  I  call 
the  approach  of  definitions  9  and  11  MAGNITUDE  SEMANTICS  and  the  alternate 
approach  SIGNED  SEMANTICS.  The  two  approaches  are  theoretically  equivalent. 
However,  since  magnitude  semantics  appears  somewhat  more  natural  and  sim¬ 
plifies  various  proofs,  it  is  the  default  for  the  rest  of  the  paper.  In  the  places 
where  signed  semantics  proves  advantageous,  it  will  be  mentioned  explicitly. 

2.3  Comparing  Two  Behaviors  over  Intervals 

It  turns  out  to  be  somewhat  more  complicated  to  compare  two  behaviors  over 
the  intervals  between  transitions.  What  does  it  mean  to  says  that  one  curve  is 
lower  than  another  over  an  interval?  To  do  pointwise  comparison,  some  notion 
of  corresponding  points  is  necessary. 

The  intuition  for  the  requisite  comparison  is  displayed  in  the  explanation  of 
spring  behavior  that  was  presented  in  section  1.1. 

If  the  mass  of  the  block  increases,  the  force  on  the  block  is  the 
same  .... 

Yet  this  doesn’t  mean  force  is  invariant  as  a  function  of  time — that  isn’t 
true.  Consider  the  time  when  the  small  block  is  at  its  rest  position;  the  spring 
applies  no  force.  Hut  since  the  large  block  is  moving  more  slowly,  it  won’t  have 
reached  the  rest  position  and  so  there  will  be  a  force  applied. 

What  the  statement  means  is  that  force  is  invariant  as  a  function  of  position. 
For  every  position  that  the  block  occupies,  force  is  equal  in  the  two  systems, 
even  though  the  two  blocks  occupy  the  positions  at  different  times.  Although 
parameters  are  defined  as  functions  of  time,  they  often  need  to  be  compared 
from  the  perspective  of  other  parameters.  Here  it  proved  advantageous  to  con¬ 
sider  force  as  a  function  of  position  Although  people  understand  arguments 
that  leave  these  changes  of  variable  implicit,  the  notion  must  be  made  precise 
and  explicit  if  computers  are  to  perform  comparative  analysis.  The  notion  of 
perspective  is  foundational 

Definition  10  A  parameter,  A',  is  railed  a  COVERING  PERSPECTIVE  over  a 
transition  interval  ( 7, ,  7, > i  )  when  the  following  three  conditions  hold: 

1.  QDlR(X,T(i,),T(nt+i))  ^  std 

2  .V||, 

»■  -ML. 

When  just  the  first  condition  holds.  X  is  called  a  partiai.  perspective 


When  a  parameter,  A  ,  is  a  partial  perspective,  it  is  strictly  monotonic  so  its 
inverse  A'-1  exists.  This  means  that  it  is  possible  to  reparameterize  any  other 
parameter,  F ,  by  composing  it  with  the  inverse: 

Fx(*)=F(  X-'(x)) 

When  A'  is  a  covering  perspective,  then  Fx  and  F\  have  the  same  domain. 
Covering  perspectives  will  prove  especially  important  in  the  inference  rules  of 
section  3. 

Definition  11  Given  a  parameter  F .  i  partial  perspective  X .  and  a  transition 
interval  (fi,  7l+i),  let  Fx  denote  F  as  a  function  of  A  .  LetU  be  the  intersection 
of  the  domains  of  Fx  and  Fx 

U  =  (A'(T(7,)),^(T(7.+i)))n(A-(f(7l)),A'(f(7l+1))) 

Define  the  RELATIVE  CHANGE  (RC)  of  F  over  ( 7l ,  7i  + 1 )  from  the  PERSPECTIVE 


of  X  as  follows: 

^flbU+l) 

ifVxeU  \Fx(x)\  >  |7*(x)| 

n  (*,.», 

if'ilEU  |F^(x)|  =  |Fa(x)| 

^(M+l) 

t/v*ec  |7a-(x)|  <  |PA(x)| 

In  other  words,  force  is  ||  from  the  perspective  of  position,  if  for  all  positions 
that  are  assumed  in  both  simulations  (Vx  E  F)  the  corresponding  forces  are 
equal.  The  definition  of  partial  perspective  sa\s  when  is  it  possible  to  use  a 
parameter  as  a  perspective;  section  3  addresses  the  question  when  is  it  useful 
to  do  so. 

2.4  Time  as  a  Perspective 

Although  comparisons  of  parameters  that  have  been  reparameterized  by  per¬ 
spectives  are  more  common,  sometimes  is  is  useful  to  compare  via  corresponding 
times.  To  keep  notation  consistent,  I  will  call  this  using  time  as  a  perspective.’ 
The  goal  is  to  come  up  with  a  meaningful  definition  for  ,+  1)  and  the  other 

RC  values. 

One  problem  is  t  hat  t  he  durat  i< ui  of  I  lie  t  wo  t  ime  intervals  might  be  different 
If  so  time  acts  as  a  partial  perspective  one  quantifies  only  over  time  in  the 
shortest  interval  Another  problem  is  that  the  two  tiansition  intervals  might 
start  at  different  tunes;  in  fact  one  interval  might  end  before  the  other  starts. 
p  %  ■  7(7,)  >  7(7,  +  ]  )  I  he  sol ut ion  is  to  align  the  inter  >.'! Is  hef< >re  quant ifv mg 


wag 


Definition  12  Given  a  parpameter  P  and  an  tnterval  (7,,  7;+i).  Lei  U  =  (0 ,  d) 
where  d  =  ML\\T(fi+i)  -  T(y, ),  T(7,+ 1 )  —  T(7,))).  Define  the  relative 
change  (RC)  of  P  over  (■>, ,  7,+i )  from  the  PERSPECTIVE  OF  TIME  as  follows: 

nL+u  Wei'  |R(T(7i)  +  f)|  >  |P(T(7,)  +  01 
P||(Ti+1)  */V<  €  U  |P(T(7.-)  +  01  =  |P(T(7.-)  +  01 
Pw£i+n  if 'it  €  U  |P(T(7i)  +  0l<|P(^(7i)  +  0l 


3  Differential  Qualitative  Analysis 

This  section  presents  ;i  numb#r  of  nil*-*  .r  ci  iiup'.it inn  atul  iii.uiipul.it mg  RC 
values,  describes  how  the  rules  wi-re  t :i-  •  >rp.  rated  into  a  computer  program, 
and  evaluates  the  program's  performance 

•  The  duration  rule  formalizes  "distance  e.pials  rate  times  time.” 

•  The  interval  derivative  rule  expresses  t  ti«-  relationship  between  one  deriva¬ 
tive  and  another  eg  m..r.  accelerat  ion  leads  to  higher  velocity" 

•  Ihe  transition  d#riv.\tiv#  rub  predicts  tie-  final  value  of  a  derivative  like 
velocity. 

•  The  self  reference  rule  says  l  hat  ever\  parameter  appears  unchanged  from 
Its  OWn  perspective 

•  The  perspec  e  flipping  rule  allows  ,  r#-a>#>ii#-r  to  change  perspectives 

•  lire  transition  and  interval  con-taut  rules  show  the  relationship  between 
constants  and  H< '  values 

•  The  end  of  time  rule  says  that  oth<  t  tilings  being  equal  a  parameter 
changes  more,  the  longer  it  is  <  hanging 

•  Idle  one  s  own  dernat tve  rule  predn  is  w hat  happens  when  a  parameter  is 
defined  m  terms  of  itself 

•  The  mult  iplicat  ion  rule  J'  nr  msl  rat<  -  t  h  c t  i  he  familiar  rules  of  qualilai  ive 
arithmetic  apply  to  HC  values  as  well  as  derivatives 

Each  of  the  rules  are  press*  tiled  as  theorems  since  they  are  proven  sound  For 
simplicity,  however,  only  the  interesting  an  I  difficult  proofs  have  been  included 
in  this  paper  Ihe  rules  have  i.e.n  implement.-, |  as  pari  of  ('A  a  7t.lAI.ISP 
program  which  solves  comparative  analysis  problems  using  |)Q  analysis  ('A 
uses  a  constraint  propagator  t>.  d-  rive  implc  iimns  ..f  t  hes*-  rules  !  In  resulting 
dependency  struct  ur.  can  I..-  t  ra  ns  l.i  t  •  -  I  into  an  I  nglish  explanation  or  us.  d 
by  an  explanation  based  gen-  r  ali/er  All  hough  (  A  is  incomplete  (there  are 
some  problems  it  for  which  it  terminates  without  sol.iugl,  it  . I o*a>  answer  and 
intuitively  explain  a  large  class  of  problems  If.  cause  the  rules  have  been  proven 
correct,  ('A  is  guarani#  is!  to  r<  .t.  b  •  .nly  sound  conclusions 

3.1  Duration  Rule 

I  Ins  rule  is  the  basis  for  tin  very  p.  m.  rt  d  ml-  pi . list  .nice  equals  rate  t  tines 

duration  If  th»-  rat#-  is  slower  in  tin  second  Simula' k.ii  then  it  will  talo  longer 
I.,  go  I  lie  sarin-  distance  Abb"  ich  ilu-  may  'i-.iii  obvious  perspectives  are 
required  to  rnak#-  pr#’#  is#-  t  li<-  ■  i  ■  •  ’  i  n  "I  rate  is  sl.wer  this  makes  il  subtle 
[i#-fore  I  ran  slat#  the  t  lie#  .|,  i . .  the  not  c  n  .  .|  distance  must  b>-  mad#  .  bar 
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Definition  13  /.ft  A  6c  a  parameter  which  ts  increasing  and  positive  (or  dt  - 
creasing  and  negative)-  over  the  transition  interval  (•*,. -)l+\).  Ihfint  DiSiAMK 
BY  A  owr  (-),■.  7j+i)  as  the  rilattve  change  of  the  distance  traveled  our  X  o i  t  t 
the  interval  as  shown  in  the  following  table : 

Starting  HC  \'alue 

ft  II  ft 


ts:. 


it 


Ending  ft 
RC  || 
Value 


?  ft  ft 

ft  II  ft 

ft  ft  y 


Note  t hat  the  parameter  A  has  a  double  purpose  m  this  theorem  it  has  1 
as  its  time  derivative,  and  it  is  also  the  perspective  from  which  V  is  seen  to  v 
In  the  following,  it  may  be  helpful  to  think  off  as  velocity,  and  A  as  position 

Proposition  1  Duration  Rule 

Lit  l  and  A  be  parami  ti  rs  such  that  A  is  a  partial  perspective  nrcr  (y, .  st  +  )  ) 
(liven  \  -  £  A.  Cft*  (+n.  and  -nis  t  anck-BY  Aft, ,  I+  , ,  the  n  T{  y1  + ,  )  -  7  { )  > 
7(  *ti+  i '  -  7(y, ).  if.  the  duration  of  y,.)|)  will  tncnasi 

Proof:  Note  that  the  proof  is  not  obvious:  \'x  f-  ^  I  prove  the  case  m 
which  tits  I  anck-BY  A||m  l+  ,  This  is  equivalent  to  requiring  A  to  be  a  covering 
perspective  Let  a  ---  A(7(y,))  and  b  -  A(7(y,+  j))  Since  A  is  a  c.oermg 
perspective.  A  lias  an  inverse  function  taking  position  to  time 

A  (a .b)  — -  (  7(y,l.  7" ( -> ,  + 1 )) 

The  function  A-1  also  exists,  has  the  same  domain,  and  a  possibly  lilbreni 
range  (7(y, )  7(y,+  1 ))  Hy  definition  V  ft/);  1+1 ,  tie  aus 

|r(A'V))|  <  |V*(  A'- ’(-r  1)1  Vi  e  (a.b) 

Consider  the  case'*  where  l  >  0.  this  implies  that  all  values  of  l  are  greater 
than  zero  because  otherwise  the  two  systems  would  have  different  transitions 
violating  the  topological  equality  assumption  fills  means  that 


(l  •  I  ( A'  <i))  <  l  (A  Vi  E  (a.b) 


S. . 


-V* 

vV 

vv 


I)  Vi  e  i  ,1  h  I 


r(A'-|(r)|  C (  A  'ir.n 

*  A  umilAr  'Himlion  iiwwl'-  f<<r  t  h#*  •  am-v  <»f  iii<  t'  .vsniK  i \  *■  aii<1  •  !■  •  i * | .. .«. » i , 

Mil*  tlrfinif  lofl  WouM  hr-  <Mftl|*W  f->  III  MKIIf'l  (.  * 

'  I  lw  .  am-  whr-re  V  <  (»  i%  Vinul.il  t  h«*ir  m  n«.  .  am-  «%  Ii«-i«  \  n  l.e,  „n«  th*  n  \  »..mH 
t*‘<v  h  a  trAmit  i*»n 


So: 


Ja  WA''l(jr))  J*  WA-'(-r) 


dr  >  0 


But  by  the  chain  rule,  the  time  derivative  of  A'  1  at  x  is  , . . —-U- -  So 

J  l  (  A  - 1 1 r  1 1 

-Y"l(6)  -  .Y-'(fi)  >  A'~  l(6)  -  A'-'(.i)  >  0 


Thus:  T( y1  +  i  )  —  ’/(*.,)  >  T(  •),  +  [ )  -  /(-,).  In  i>i  her  words,  the  d  u rat  ion  of  t  he 
interval  increases.  □ 

It  would  be  nice  if  one  could  show  that  the  duration  rule  was  sound  if  the 
premise  was  weakened  to  have  1+ll  for  some  arbitrary  covering  perspective 
P.  However,  the  following  proposition  shows  that  this  is  false;  just  because 
P||*i+1-  for  a  perspective  .Y  doesn't  mean  that  there  doesn't  exist  some  other 

perspective  Z  such  that  Pftf,  I+ll 

Proposition  2  Noli  Uniqueness 

Given  a  system  u  ith  paramete  rs  P.X.Y.  and  Z  such  that  A  V  and  Z  art  cov¬ 
ering  perspectives  over  (y,.  -,1+1).  then  it  is  possible  that  Pft*  1+|,  and  P\\'t  l+)) 
and  P U, ,  ,+  |  | 


The  example  shown  in  figure  7  ill  list  rates  t  lie  proof  bv  construction  The  thin 
lines  indicate  the  values  of  the  first  system  while  the  dotted  lines  indicate  the 
value  of  the  second  system  I  lie  first  row  shows  that  from  the  time  perspective 
the  behavior  of  P  doesn’t  change  1  he  second  row  shows  the  relative  change  of 
the  perspectives  lie  third  row  depicts  P\  /’>  and  f’z 

Alt  hough  t  his  as  pee  t  of  H<  values  may  seem  st  range,  it  is  act  ually  inevitable 
After  all.  everything  is  relative  tootles  perspective  Imagine  a  machine  which 
hourly  logs  the  linearly  increasing  concent  rat  ion  of  alcohol  in  a  fermentation 

tank  It  produces  the  following  sc, |uciu  e  of  measurements  0  02  (I  04.  0  0(>.  0  (IK. 

etc  Hut  m  the  identical  tank  m-arby  tie  |.  ggmg  ma>  Inn,  has  a  defective  motor 
which  runs  loo  slow  ]y  and  delay  s  i  lie  u  |,  asur  <  i  neiil  s  A  It  hough  the  fern  ten  tat  toil 
is  proceeding  at  the  same  pace  m  both  tanks,  the  second  log  will  read  0  03. 
0  06  0  09.  0  13  etc  1  luis  the  plant  inspector  win,  only  sees  the  alcohol  tune 
urve  from  t  lie  persp.ctu,  .|  th*  logging  T  vie  might  think  that  second  tank 
was  fermenting  nmo  pm  kly  even  though  tin-  only  r«  a  I  change  was  a  slowdown 
in  the  speed  of  the  timing  im  lor 


3.2  Derivative  Rules 

I  lies,  rules  .  onno,  i  par  am- tore  t  h  0  if  iii:„  Privative-.  |  he  first  works  over 
intervals  and  *  In*  second  predicts  IP  values  it  mtei v.al  endpoints  I  he  tut  nit  ion 
1 1 1  hind  the  ti  rst  I-,  if  ,i  |  ■  ,  r  a  m»  i  *  i  i-  at  1 1,*  -tot  if  an  interval,  but  itsderivative 
|s  .>  1  ver  t  In  interval  1 1,>  n  1 1,,  par  , m>  o  i  . -I  1  ,  ,v  >  r  l  In  inter va I  As  always 


i  ' 


Figure  7:  Pft(A0  ,,  yetP\\'oi)  andP^Q  ,, 


the  ubiquity  of  perspectives  complicates  the  matter  Note  the  special  role  of  X 
both  as  perspective  and  second  integrand  of  .4 

Proposition  3  Interval  Derivative  Rule 

Let  .4.  V',  and  X  be  parameters  such  that  A  =  V',  V  =  Jj.Y.  and  X  is  a 

covering  perspective  over  (7e7i  +  i)  Furthermore  let  .-1  and  V  be  positive  ,ner 
the  interval  (7,,  7,+  1 )  //-d'ft,,  .411,.  and  .4^*,+  1)  Then  I'U* , 

Proof:  The  chain  rule  makes  this  rule  considerably  harder  to  prove  than  the 
duration  rule  It  suffices  to  show  that  there  exists  some  position  such  that 
|  <  1  for  all  positions  up  to  and  including  this  position  Once  it  is  known 
that  V  goes  down,  the  same  argument  can  be  used  to  show  that  it  continues  to 
go  down  Thus  it  will  stay  down  until  7,  +  1  is  reached 
I.et 


dt 

'u,  =  d7  = 


V(.Y->(*)) 


Hx)  _  £i 

dz2 

.4  can  be  expressed  as  a  function  of  .Y 
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•  If  t he  object  is  accelerating  slower,  then  its  terminal  velocity  will  he 
smaller. 

•  If  the  object  is  accelerating  at  the  same  rate  hut  starts  with  a  slower  initial 
velocity,  its  terminal  velocity  will  be  smaller 

•  If  the  object  accelerates  for  a  shorter  distance,  then  it  will  finish  going 
slower. 

Proposition  4  Transition  Derivative  Rules 

Let  .4.  1‘.  and  X  be  parameters  such  that  .1  =  j.  V  —  jL  \\  \  ,s  a  partial 
perspective  over  (■),.  -,1+1 )  and  both  .4  and  \  are  positive  over  the  interval.  If 
one  of  the  following  conditions  is  time. 

•  Vlf,  and  .4||*  l+1,  and  DISTANCE-BV  .Y||( ,  1+n 

•  (,.41j1  A  .4U*,+  n;  and  I'U,  and  D1STANCK-BV.\'||(I  l  +  n 

•  DISTANCE-BV.\ l+ 1 1  and  Iff,  and  .4||*1+1) 
then  n;i+1. 

The  rule  is  quite  a  mouthful,  but  that  is  simply  because  it  is  very  general. 

3.3  Perspective  Rules 

These  rules  deal  with  establishing  RC  values  for  perspectives  and  switching 
between  them.  The  first  is  very  simple,  but  turns  out  to  be  quite  important 
The  intuition  is  that  if  the  plant  manager  was  foolish  enough  to  try  and  use  the 
logging  devices  to  log  their  own  speed,  he  wouldn’t  get  a  useful  result  Both 
the  normal  and  slow  machines  would  record  that  they  turned  one  full  revolution 
during  each  revolution  of  the  timing  motor. 

Proposition  5  Self  Reference  Rule 

For  any  parameter  P,  if  P  is  a  covering  perspective  ore  r  (■),,  ->1  +  1)  then  /’||(|1+11. 

lhe  perspective  flipping  rules  switches  between  perspectives.  The  intuition 
is  that  flipping  perspectives  (i.e.,  X  f'  to  Px  )  flips  ft  to  fi-  if  both  parameters  are 
positive  and  increasing  over  the  interval. 

Proposition  6  Perspective  Flipping  Rule 

If  the  parameters  X  and  P  are  valid  perspective  s  ore  r .  *.  1  + 1  ).  the  sign  of  A 
equals  the  sign  of  I'  over  the  interval,  and  A  ft  ^  I+Il.  then 

/’ft,',*,,  ifQDIH(. Y.7(-„).7b,+i))  /  QW/f(/>.7<y,).7(*,+  ,)) 
Hl*>+I)  »/<^/W(.Y.7b,).7b1  +  1))  =  /*.  7b,  ».7<  •>,♦,!) 

If  the  sign  of  X  is  the  opposite  of  the  sign  of  P  then  the  RC  values  are  reversed 


Proof:  I  will  prove  the  ca.->e  where  hotli  V  and  P  are  increasing;  the  other  cases 
are  almost  identical.  Let  a  =  .V  (T(->,)).  and  b  =  A’(T(7,+  1)).  For  an  arbitrary 
■r  €  [a,  b)  3p  such  that  A  (  P~  1  (p))  —  x  because  P  is  a  covering  perspect  ive.  and 
thus  onto.  Let  t i  =  P~  '(/<),  and  let 

?=  A'«.)  =  A  t  P~ 1  (/') 

By  the  definition  of  A  ft^  1+ ,  ,  it  follows  that  x  >  x.  Let  to  —  A’_1(J‘)-  Since 
A’  is  increasing  to  <  P  Again  because  P  is  onto.  3p  such  that  P~x{f)  =  / 0  so 
X{P~x(p))  =  x.  Now.  f  <  p  because 

P~\F)  =  to  <  >i  = 

and  P  is  increasing  But  this  means  that 

P(.Y~V))  <  P(X~'(x)) 

and  since  r  was  arbitrary,  it  follows  that  /’ij,'  1+|)  □ 

3.4  Constants 

Frequently  a  system  will  contain  a  few  constant  parameters  whose  values  never 
change.  The  following  rules  are  a  simple  way  to  express  relationships  between 
constants  in  the  notation  of  comparative  analysis.  The  intuition  is  that  since 
perspectives  just  scale  time,  and  constants  don’t  change  over  time,  all  perspec¬ 
tives  agree  on  the  behavior  of  constants  If  there  was  no  fermentation  happening 
in  either  vat  (i.e.  the  alcohol  concentration  was  constant  in  both  vats),  and  the 
concentration  of  alcohol  was  higher  in  vat  two,  then  both  logging  devices  would 
agree  on  this  even  though  their  timing  motors  differed. 

Proposition  7  Transition  Constant  Rule 

If  a  parameter  K  is  a  constant  over  (■>,.  -ft  +  1 ),  and  A’ ft;  then  A'fti+1. 

Proposition  8  Interval  Constant  Rule 

If  a  parameter  K  is  a  constant  over  (7, ,  y1  +  1 ).  and  A' ft,  then  for  all  parameters 
P,  if  P  is  a  covering  perspective  over  the  interval  (7v,7,>i).  then  K  ft(^  1+1). 

3.5  Rules  with  Time  as  a  Perspective 

It  is  very  common  for  one  parameter  to  be  the  derivative  of  another  with  respect 
t  o  t  ime  When  it  is  possible  to  reason  abi  >ut  t  hese  relations  from  t  he  perspec  I  ive 
of  time,  greater  power  is  achieved  because  the  chain  rule  doesn’t  interfere  as  it 
does  in  the  derivative  rule.  The  only  drawback  is  the  fact  that  tlu'se  rules  are 
less  frequently  applicable. 

The  first  rule  says  that  if  the  a  parameter  is  j[ fr< >in  the  perspective  of  time, 
and  the  duration  of  the  interval  is  increasing,  then  the  parameter  will  have 
changed  mom  by  the  etui  of  the  niter,  il 


p  i  ■  •  j  — it  it  >u  •*  rh.  k  11.1  i>i  i  mi*  n  nit- 

'  \  ’■  ;  ;  fi  I  \  _  <1  n  •/  A  \lt,-  /•'  /  A  ftf  Iht  >(■///  <</  A 

■  •  -  •  ■  J  ■  . :  'I  r,  |.|  ■  •  li  !  *• .  -■ ,  ..  ■  i  ,1 II  .1  .1  hi  th(  slt)H  i  f  X  <  d<  lllntll  l  If  I  III 

:  .  .  .  '  •  .  - . .  i>  ••  iX  i, 

\  i '  i  -  ■  1  ‘t '  i'ii  i 

'  "  j  i  - 

:  ’  >  ?  I 

",  :  r  l  t  t  !,i>  >•  mu n.t  i'  iruml  and  t  litis  '  unit  i 1 .  i ■  u I  it  should  be  not 1 

i  ■  t-.- r  t  •  \|  r>  ->  using  sinned  semantics  I  In-  second  rule  is  used 

"  •  r : r ^  t  paratii'  f ♦- r  H<‘  value  fr-  >m  the  perspective  r  f  time  It  applies 

a  :■  tiin<  •  )»  m  itive  >f  a  parameter  is  a  linear  function  of  tin*  parameter. 

I ’ rt >(•> >'it i. >n  10  One's  Own  Derivative  Rule 

■  1  '  .  ;  :  i  l\  h  |•■|| 'ini  I' I"  'Ui  h  Ihul  \  —  -77  \  .  1  =  Ml  I,  I(  A  .  /\  ) .  illlilh 

.  -  1  I,  •  j  1  rr  ■  .  ml. nil  If  \  ■  1 1  *,  >  1  S  1 1  ,1  tl  d  l  (lit  il  l\  ,  1+  ,  that  +  ,  , . 

3.R  Halt's  from  Qualitative  Arithmetic 

s.  1; ,  1  i,,t.ti  n-  m!  01.  tii  7  L't.  1  V  li.ts  <|e veil. ped  constraints  on  deriva- 
•  ■  ■  v.«'.n«*s  !  r  parameters  m  M0>  mi  1. 1  and  in.  uotouic  function  constraints. 

I  -  -  \  tt-.pe  0  \  •  1  -  iti- 1  t li**  derivatives  <>f  A  and  )’  are  positive,  then  7. 

.«•  1  if  ;  stove  i.  r  1  v  it  1  \ .  a>  v\e||  fliese  rules  can  he  generalized  to  include 

h  ■  .  »  ,>-s  »t  transit).  11  |». ints  an-l  ever  intervals  II >  re.  I  present  just  the  rule 

i  M  !  list  ran .  1  it  a  1  r  in-it  1.  11  point 

Proposition  11  M  ult  iplii.at  ion  Rule 

'i  i>,;.’  1  n  ;.'i'.i'«.'.>>  nhoh  ’in  rthih ./  hi/  tin  ron.it  rant  1 ,  7  ~  Mt'I.T(A.  V). 
.  '  n)  lilt’  Ini  If/-  tin  />.  -  >  >  1  /•  /.  /fi  '  1, linn  fur  7  nl  ,1  transition 
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3.7  Implementation 

To  test  the  theory  of  HQ  analysis,  a  program  called  CA  has  been  written  on  a 
Symbolics  lisp  machine.  When  a  user  selects  an  example,  (TV  runs  QS1M  [15] 
on  the  example  to  produce  a  set  of  qualitative  behaviors  for  the  example.  The 
user  selects  a  behavior  and  also  a  set  of  initial  RC’  perturbations.  C'A  translates 
the  QSIM  behavior  and  perturbations  into  ARK-4  assertions.  At  this  point  ARK 
forward  chains  using  the  propositions  described  earlier  in  this  section. 

Each  of  these  propositions  is  implemented  as  an  ARK  rule  or  more  than 
one  if  the  proposition  used  disjunction  or  negation.  For  example,  the  duration 
rule  (proposition  1)  is  encoded  as  the  three  ARK  rules  of  figure  8.  The  various 
definitions  and  propositions  require  about  sixty  ARK  rules. 

( =>  (AND  (d/dt  ?x  ?v)  ;  ?v  is  the  derivative  of  ?x 

(distance-by  ?x  (?start  ?ettd)  deq)  ;  ?x  travels  the  same  distance 
(RC  ?v  ("’start  ?end)  ?c  (P-  °x))  :  the  RC  of  ?v  is  ?c  from 

(opposite-RC  ?c  ?oc))  ;  the  partial  persp.  of  °x 

(duration  (?start  ?end)  :'oc)  ;  if  ?c  is  ff,  ?oc  is  JJ. 

duration-rulel ) 

(=>  (and  (d/dt  ?X  ?v) 

(distance-by  ?x  (?start  ?end)  °oc)  ;  if  ?x  travels  ?oc  distance 
(RC  ?v  ("’start  ?end)  ?c  (P-  ?x))  ;  and  V’s  RC  agrees 

(opposite-rc  ?c  ?oc)) 

(duration  (?start  ?end)  ?oc)  ;  then  the  duration  is  ?oc 

duration-rule'2) 

(=>  (and  (d/dt  x  ?v) 

(distance-by  ?x  (?start  ?end)  ?oc)  ;  if  ? x  travels  ‘less"  distance 
(RC  ?v  (?start  '’end)  deq  (P-  ?x)))  ;  and  V  doesn’t  change 

(duration  ( ? start  ?end)  ".’or)  ;  then  the  duration  is  ‘less’ 

duration-ruled) 

Figure  8:  Propositions  Are  Encoded  Directly  Into  ARK  Rules 

The  simplicity  of  the  transformation  from  proposition  to  ARK  code  provides 
confidence  in  the  soundness  of  the  implementation.  And  the  fact  that  most  rules 
get  used  in  each  explanation,  establishes  their  utility. 

Since  it  tsan  initial  prototype,  ("A  makes  no  use  of  control  rules.  All  possible 
forward  chaining  inferences  are  made  using  every  possible  perspective.  Despite 
this,  computation  rarely  exceeds  a  minute  on  anv  (1f  t|1P  problems  tested.  If 

4  A  R  K  is  a  (Ifsccti'Iiiiil  .  >  f  A  MOB  I  >  ['.  j  inij>T‘ni.  nleT  !,v  Itnwir  Slmil.r  *md  oilier* 


larger  problems  were  to  be  attempted,  some  form  of  control  would  be  desirable. 
Backward  chaining  from  a  goal  pattern  might  increase  efficiency.  There  appears 
to  be  no  reason  why  the  schemes  of  [4,2-1]  could  not  be  applied.  Possible  heuris¬ 
tics  include  preferential  investigation  of  certain  perspectives  and  avoidance  of 
certain  computationally  explosive  rules  like  the  perspective-flip  rule. 

Another  technique  to  speed  up  reasoning  is  explanation  based  generalization 
[16,5].  Following  the  approach  of  [13],  I  implemented  a  postprocessing  learning 
routine  that  takes  CA  explanations  and  produces  new  ARK  rules  which  may  be 
added  to  the  ones  presented  above.  While  these  new  rules  are  independent  of 
any  particular  domain  (i.e.,  springs),  they  are  optimized  to  solve  a  specific  class 
of  comparative  analysis  problems.  Less  general  than  the  rules  presented  above, 
the  new  rules  are  considerably  more  general  than  the  specific  explanation  from 
which  they  are  derived.  Although  I  have  completed  the  FBG  implementation, 
the  empirical  evaluation  of  EBG's  ability  to  increase  DQ  processing  efficiency 
remains  as  an  area  for  future  research. 


3.8  Differential  Analysis  Suffices  for  Most  Examples 

Since  ARK  maintains  justifications  for  all  its  assertions,  it  is  possible  to  generate 
explanations  for  CA’s  conclusions.  Consider  the  spring/block  example.  The 
question  here  is:  “What  happens  to  the  period  of  spring  oscillation  if  the  mass 
of  the  block  is  increased?”  The  system  is  defined  in  terms  of  six  parameters: 
spring  constant  — K ,  mass  M,  position  X,  velocity  V,  acceleration  A,  and  force 
F  obeying  the  following  equations: 


F  =  mult(A/,  A) 

F  =  mult(— A’,  A') 

±M  =  sid 

21  -I<  = 

The  initial  conditions  are  specified  as  follows:  A/(0)  >  0.  — /\'(0)  <  0,  V ' (0 )  = 
0,  and  A’(0)  =  x0  <  0.  Since  energy  conservation  is  not  made  explicit  in  the 
equations,  QSIM  produces  several  possible  behaviors  for  this  system.  Although 
comparative  analysis  could  be  done  on  any  of  the  behaviors.  1  assume  in  this 
example  that  the  user  selects  the  interpretation  corresponding  to  stable  oscilla¬ 
tion. 

Now  the  user  selects  the  perturbation  Because  some  parameters  depend  on 
one  another,  not  all  parameters  may  be  perturbed  The  situation  is  analogous 
to  the  problem  of  specifying  a  unique  solution  to  a  differential  equation  where 
values  must  be  given  for  the  independent  parameters  and  a  set  of  boundary 
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conditions  provided  In  t Ins  example.  M  and  -K  are  independent,  while  valm- 
for  A  and  l’  are  needed  as  boundary  conditions  1'hus  to  specify  a  comparative 
analysis  problem,  these  four  parameters  u<  ■  d  to  he  given  initial  RC  values  ' 
For  this  example,  the  perturbation  consists  of  the  following  initial  RC  values' 

m0-  -^Ito;  Vllo-  a,1(1  A'H." 

Given  this  input,  CA  correctly  deduces  that  the  block  will  take  longer  to 
reach  the  rest  position  (.V  =  0)  from  its  original  negative  stretched  position. 
Figure  9  shows  the  explanation  that  CA  generates;  this  is  created  by  throwing 
away  alt  perspective  information  om  e  computation  is  finished.  1  have  annotated 
the  explanation  with  tin  names  of  rules  used  in  each  step. 

Assuming  M  is  increased: 

X  doesn't  change  and  ( self-reference  rule) 

K  doesn’t  change  and  (interval  constant  rule) 

F  equals  -K  times  X 

So  F  doesn’t  change  ; 'multiplication  rule) 

and 

M  increases  and  (interval  constant  rule) 

F  equals  M  times  A 

So  A  decreases  (multiplication  rule  } 

So  V  decreases.  (derivative  rule) 

So  the  time  duration  increases,  (duration  rule ) 

Figure  9;  ('A  Generated  Kxplanation  for  Spring  with  Heavier  Block 

At  present  CA  has  been  tested  for  multiple  perturbations  on  over  a  dozen 
examples  including  the  RC  circuit  shown  in  figure  10  While  if  always  terminates 
and  never  produces  an  incorrect  answer.  CA  doesn’t  necessarily  deduce  R(  ’ 
values  for  every  parameter. 

3.9  Differential  Analysis  is  Incomplete 

As  is  explained  in  the  sections  below,  different  types  of  ambiguity  are  the  cause 
for  the  incompleteness  of  HQ  analysis  Vet  J)Q  analysis  handles  ambiguity  dif¬ 
ferently  from  other  forms  of  qualitative  reasoning  For  example,  when  QSIM  is 
faced  with  ambiguity  about  a  parameter’s  value,  it  branches,  spawning  perhaps 
three  new  behaviors:  one  with  the  parameter  equal  to  a  landmark  value,  one 

"'The  choice  of  these  four  parameters  is  somewhat  arbitrary  Matheinat ieally.  it  would  he 
equally  reasonable  to  choose  .-l  instead  of  .V,  but  this  does  not  make  physical  sense;  it  seems 
intuitively  impossible  to  directly  affect  acceleration.  Since  there  is  no  wav  to  deduce  ibis 
from  the  differential  equation  model,  it  is  esenfial  for  the  person  who  constructs  the  model  to 
annotate  the  structural  description  with  the  list  of  4 causally  primitive’  parameters  in  this 
case,  the  four  listed  above. 


Figure  10:  An  RC  Circuit  Based  on  the  Wheatstone  Bridge 

greater  and  one  less.  QSIM  can  do  this  because  the  nature  of  inequality  guar¬ 
antees  that  either  A  <  B,  A  —  B,  or  A  >  B  While  this  is  true  for  RC  values 
at  transition  points,  it  is  not  true  for  RC  values  over  intervals 

Proposition  12  Non- Exhaustivity 

(liven  two  parameters,  V'  and  P,  such  that  P  is  a  covering  perspective  over  an 
interval  (7,.  7,+  i  ),  it  is  not  necessarily  the  case  that  one  of  1+1,.  1||(^ 
or  l  .+  n  holds 

Proof:  Appendix  A  provides  an  example  which  proves  this  statement  □ 

Thus  unlike  qualitative  simulators,  DQ  analysis  can  not  branch  when  faced 
with  uncertainty,  it  simply  acts  mute  The  following  sections  explain  the  three 
lactors  that  can  cause  DQ  analysis  to  fail  to  predict  all  of  the  relative  changes  in 
a  perturbed  system:  ambiguous  questions,  ambiguity  resulting  from  the  quali 
tative  arithmetic,  and  the  lack  of  a  useful  perspective 

3.9.1  Ambiguous  Questions 

Some  questions  simply  don’t  contain  enough  information  For  example:  "What 
would  happen  to  the  period  of  oscillation,  if  the  mass  of  the  block  was  heavier 
and  the  spring  was  more  stiff?”  There  is  no  answer  to  this  question  because 
it  is  inherently  ambiguous.  The  increased  mass  tends  to  increase  the  period, 
but  the  increased  spring  constant  tends  to  decrease  it  Thus  the  duration  might 
increase,  decrease  or  remain  unchanged 


3.9.2  Ambiguity  Introduceil  by  Qualitative  A  i  it  lintel  ic 

Since  F>Q  a  Italy  sis  uses  the  sane-  -in  i  ,ii.t  i  iv  >  ant  luu<  1 1>  nl  ih/ed  by  ■  »t  b*  •  r  t>  this 
of qualitative  rea-onum  it  Meadd  nd  F.-  surpriMiig  that  ambiguity  r.uw's  a 
problem  here  as  well  I  In  pr>  •Mein  is  f  ..  a.  ,1  m  lie  lac  i  i  hat  qtialit  at  iv values 
p>f  winch  R( '  values  are  an  instance)  do  icl  f  irm  a  group  over  addition  [lb] 
As  a  result,  unique  inverses  do  not  always  .  mMs  ami  it  is  frequently  impossible 
to  determine  the  qualitative  value  of  a  parameter 

For  example  consider  i  In  -  sprinn  I  loch  system  ,  ,f  t  ne  last  sec i  ion  HQ  analy¬ 
sis  r or rec tlv  predicts  that  i In-  I >|. ..  k  will  1  a k>  longer  to  reach  t In-  first  transit  ion. 
the  block’s  rest  position  But  t In-  period  of  os.- ilia' mu  requires  four  transitions 
starting  from  a  negative  initial  position.  V  moves  to  zero,  then  to  a  positive 
maximum,  then  to  zero,  and  linally  to  its  original  position  Because  ot  ambigu¬ 
ity  in  the  extreme  positions  of  .V.  I)Q  analysis  can  make  no  prediction  about 
duration  of  these  last  three  transition  intervals  Why  is  this?  Because  of  the 
qualitative  arithmetic,  it  is  impossible  to  show  that  .V |!  i  e..  that  .V  sweeps  out 
the  same  distance  when  the  mass  is  increased  Because  of  this.  X  is  not  known 
to  he  a  covering  perspective  so  the  derivative  and  duration  theorems  can  not  he 
used.  Thus  there  is  no  way  to  determine  the  ]{('  value  for  the  whole  period 
This  problem  is  directly  analogous  to  QSIM's  prediction  of  spurious  behav¬ 
iors  [15].  Given  a  Hooke's  law  description  of  the  spring/block.  QSIM  produces 
many  possible  behaviors  in  addition  to  the  correct  description  of  stable  oscil¬ 
lation.  Furthermore,  the  [)Q  problem  can  be  alleviated  in  the  same  way  that 
kuipers  caused  QSIM  to  disregard  behaviors  other  than  stable  oscillation-  bv 
augmenting  ihe  structural  description  with  equations  describing  conservation 
of  energy  Now  GA  can  deduce  that  since  potential  energy  is  equal  to  force 
times  distance,  increasing  the  block's  mass  leaves  total  energy  unchanged.  This 
a'lovvs  it  to  recognize  .V  as  a  covering  perspective  and  deduce  that  the  duration 
increases  for  each  of  the  period's  four  transition  intervals 

3.9.3  No  Useful  Perspective 

Other  qt.  -st ions  are  even  more  dilTicult  to  answer  "What  would  happen  to  the 
period  of  oscillation  if  the  initial  displacement  is  increased7"  Since  people  have 
trouble  with  this  question,  it  should  not  be  surprising  that  DQ  analysis  cannot 
answer  t lie  question  either.  In  fact,  the  answer  is  "period  does  not  change",  but 
the  only  way  t  o  show  t  his  is  to  solve  t  lie  d  life  rent  id  equal  ion  f  >r  an  equat  ion  for 
period  and  notice  that  it  is  independent  of  amplitude  I  he  difficulty  is  rooted 
in  the  fact  that  no  nstful  prrspn  tm  mst\  to  provide  a  handle  on  the  problem 
1  here  is  no  system  parameter  /’  Meh  that  \'1|(0  n  f'learlv  X  won't  work  as  a 
perspective,  since  it  doesn’t  sweep  out  the  same  range  m  the  two  cases  In  fact, 
it  is  easy  to  prove  t  hat  no  artificial  perspect  jve  could  satisfv  the  equation 


Proposition  13  Given  the  definition  of  l'  as  specified  above  for  the  spring /block 
example  with  .V ff 0 ,  /ct  =  7(7i)  and  t\  =  T(7i)  There  are  no  continuous, 
real  valued,  functions  P,  P  such  that 

P( 0)  =  P( 0)  =  f,()  A 
P(M  =  P(t\)  =  V\  A 
V(p-'(v))  =  V(P-i(p))  Vp€  (po-Pi) 

Proof:  Since  A’ff0 .  initial  potential  energy  is  higher  in  the  perturbed  system, 
so  kinetic  energy  is  greater  at  7^  This  means  that  V  i  e.. 

V  (<i)  =  t'i  >  t’i  =  V  (/,) 

Because  V  and  V'  are  continuous 
lim  l '(f)  =  vj 

(—(7 

and 


Pr- 

[• 


lim  V'(f)  =  t’i 

i-i, 

Similarly, 

lim  P-’(p)  =  fh 

P^Pl 

and 

lim  P'Hr)  =  Pi 

p—p  1 

Thus 

lim  (V(P~'l(p))  -  V(P~1(p)))  =  vi  -  vi  /  0 

p—pi 

So  there  exists  some  q  €  (po.Pi)  such  that 

V'(P-‘(«))  =  V(P^(q)) 

□ 

Thus  there  is  no  function,  P ,  that  can  act  as  a  perspective  such  that  V'||(0  )r 
This  really  shouldn’t  be  very  surprising.  After  all,  the  block  really  does  move 
faster.  The  only  reason  that  the  period  is  unchanged  is  that  the  increased 
velocity  is  exactly  counterbalanced  by  the  increased  distance  the  block  must 
travel.  It  would  be  foolish  to  try  and  claim  the  velocity  doesn't  increase  when 
it  does.  Instead,  an  intuitive  explanation  should  account  for  the  balance  of  the 
change  in  velocity  and  distance.  This  type  of  explanation  is  outside  the  realm 
of  DQ  analysis,  and  probably  beyond  the  abilities  of  qualitative  physics  as  well. 


3.10  Extensions  for  Diagnosis 

A  natural  application  for  coinparat  o '■  analysis  is  1 1 1*  •  automated  diagnosis  of 
continuous  devices  such  as  .mil.  ,g  ci  r-  me  or.  mts  ( >>  aerate  and  test  i>  a 
standard  paradigm  for  hardware  diagnosis  ;i'  candid. it<-  faults  ar<  proposed 
then  evaluated  to  sc-  if  tln-v  account  for  tin-  fault)  measurements  I. ike  all 
forms  of  comparative  analysis.  IX j  analysis  •  an  he  us.  d  to  test  any  candidate 
faults  that  can  he  described  as  perturbations  of  .  out muons  parameters  in  the 
device.  In  addition,  however.  DQ  analysis  has  tin  p. »t«-nt lal  to  generate  classes 
of  candidate  faults.  I  In-  key  is  to  run  tin  Hi  2  inference  rules  in  reverse  lor 
example,  the  duration  rule  says 
Proposition  1  Duration  Rule 

I.et  V  and  X  be  parameter >  sui  k  that  X  a  partial  pirf.pt  i tire  ore  r  ( 1 1  | 
Cinen  V’  =  ji  X .  and  -UMAVci  nv.Vi]11I+ll  then  the  duration  of 

(fo'fi  +  l)  til//  mere, iff 

A  natural  question  is  "Can  the  duration  rule  be  reversed'.’  Is  the  converse 
sound'7” 

Conjecture  14  Converse  Duration  Rule 

Let  V  and  X  he  para  me  tf  rs  (been  1  =  ^  .V  and  — ■  I  >  I  s  I  A  N  CK  HY.YfJ(11+ll  If 
the  duration  of  (•>,.  7,+  ;)  If .  thin  l'ff*1+1  . 

Unfortunately,  the  converse  is  false,  as  are  tin  converses  for  other  important 
rules  such  as  the  various  derivative  rules  The  problem  results  from  an  implicit 
closed  world  assumption  used  in  reversing  the  rule  that  one  of  the  three  R( 1 
values,  ff,  ij.  or  ||.  always  applies  Proposition  12  showed  that  this  was  false. 

Of  course  every  transition  interval  could  he  broken  into  pieces  such  that  a 
single  RC  value  applies  over  each  piece,  hut  this  misses  the  fundamental  issijp. 
The  decomposition  of  time  into  transition  intervals  is  forced  by  behavior  of  the 
system.  Thus  transition  intervals  have  genuine  ipialit.it ive  importance.  While 
sometimes  useful,  decomposing  transition  intervals  into  smaller  pieces  runs  the 
risk  of  introducing  irrelevant  distinctions. 

Although  the  converse  of  the  duration  rub-  is  not  sound,  its  converse  might 
still  be  profitably  used  as  a  heuristic  candidate  generator.  By  reversing  the  DQ 
inference  rules,  it  may  be  possible  to  provide  focus  to  the  search  for  probably 
faults  in  misbehaving  analog  circuitry 


4  Changes  in  Behavioral  Topology 

Hivall  lh.it  tin'  inliritiii-  ruii>  nl  yrtimi  li  r»- 1 1*  -<  1  <ui  tin  ,i»umpt  t- -n  t  li.it  t !.( 

|  '•  rt  11  r  l"'il  In'lia  Vi.  >r  w  .i>  t.  >| . .  -i.  mi.  alls  j  nal  t  defimt  i.  ■  1 1  •»  i  t  • .  tin'  mil  m  I  l.<  I,  i  \  .  • .  i 
lu  I'lln-r  words,  it  was  assumed  that  while  tin  |n  rt urliat  i<  -n  ttiiiilit  change  1 1 1  < • 
ri'lat  l  Vi-  values  of  parameters  and  st  ret  eh  or  shrink  tin-  h-ngt  h  of  t  nne  interval' 
the  underlying  sequence  of  transitions  would  not  change  Yet  pert  i j r I ■ . 1 1 n  .i.s 
often  will  change  the  order  or  nature  of  transitions  1  his  sort nm  explains  how 
to  recognize  the  changes  and  predict  the  resulting  In-havior 

lo  illustrate  these  computations.  I  use  the  simple  example  of  the  le  at  •  ■  x 
changer  (figure  -I)  from  section  1  2  1  Ins  system  is  descrihed  in  terms  of  tiv< 

parameters,  each  a  function  of  time  heat  Q.  heat  flow  F  thermal  conductivity 
-A-,  velocity  of  the  liquid  through  the  pipe  1‘.  and  position  of  a  unit  v..|um>  ■.  t 
oil”  A  The  following  equations  are  obeyed  ' 

v=  &  * 
r  =  rtQ 

F  -  mi-lt(Q.  - A)  (H 

In  addition  V'  and  —A  are  considered  independent  and  assumed  constant 
over  time.  The  initial  conditions  specify  the  value  for  tin'  independent  parame¬ 
ters:  lr(0)  >  0  and  —  A(0)  <  0.  and  also  the  boundary  conditions:  A"(0)  -  r(,  < 
0  and  Q(0)  =  170  >  0  From  this  information  the  initial  value  of  the  dependent 
parameter.  F.  can  be  determined,  denote  F(0)  —  /0  An  invariant  specifying 
that  A”  must  always  be  less  than  or  equal  to  zero  ends  the  simulation  when  the 
liquid  individual  leaves  the  pipe. 

Given  this  description,  QSIM  (and  other  qualitative  simulators  [6])  produces 
the  tree  of  qualitative  states  (state  trf.f,)  shown  in  figure  11.  Since  each  path 
through  the  tree  is  a  topologically  distinct  behavior,  this  tree  represents  three 
possible  behaviors  for  the  heat  exchanger.  The  topmost  path  (QS1,  QS2.  QSH) 
corresponds  to  the  behavior  of  figure  12  in  which  the  system  reaches  thermal 
equilibrium  just,  as  the  oil  leaves  the  exchanger. 

Because  of  its  qualitative  representations,  QSIM  cannot  choose  between  the 
different  behaviors  for  the  heat  exchanger:  as  far  as  QSIM  is  concerned,  they  are 
all  plausible.  Since  DQ  analysis  works  relative  to  a  single  behavior,  one  path 
through  the  tree  must  be  chosen  before  running  the  rules  of  section  II  This 
selection  of  a  behavior  is  a  modeling  decision;  I  assume  that  it  i.s  done  by  a 
human  The  selection  consists  of  a  series  of  choices  at  each  branch  in  the  tree. 
By  ruling  out  possible  behaviors,  each  choice  implicitly  constrains  the  model  of 
the  system,  restricting  the  possible  real  values  associated  with  the  qualitative 
values  of  each  parameter  Thus  the  selection  of  behavior  (QS1.  QS2.  QSH) 

For  simplicity,  the  simplistic  ‘liquid-individual’  model  t>f  fluids  is  used  here;  see  [to]  f.„  * 
discussion  of  the  problems  with  this  model. 

7  For  simplicity,  this  model  does  not  distinguish  between  temperature  and  heat. 
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Figure  11  QSIM  State  Tree  Generates  Possible  Behaviors 

makes  implicit  assumptions  about  the  relative  values  of  fluid  velocity.  V'.  and 
thermal  conductivity,  -K 

These  implicit  constraints  are  equivalent  to  the  unambiguous  selection  of  the 
initial  behavior  However,  the  comparative  analysis  perturbation  can  weaken 
the  balance  of  constraint  in  two  ways 

•  The  initial  behavior  can  be  rendered  inconsistent  Section  4.1  explains 
how  the  conflict  is  recognized  anil  a  new  consistent  path  is  found 

•  Alternate  behaviors  may  become  consistent  Section  1  2  explains  how  to 
locate  other  consistent  paths  through  lie-  state  tree 

4.1  Initial  Behavior  Inconsistent 

Suppose  someone  selected  the  path  (QSl.  QS2  QS3t  as  the  heat  exchangers 
initial  behavior  (figure  12)  and  chose  the  perturbation  -A" 'O'  The  state  QS3 
dictates  the  two  transitions.  Q  reaching  zero  and  V  reaching  zero,  in  the  same 
time  instant  Since  the  perturbation  causes  heat  to  be  lost  more  rapidly,  QS3 
can’t  be  part  of  the  final  behavior  If  one  assumes  that  it  is.  the  duration  rule 
(section  3  1)  deduces  a  contradiction,  aw  follows 

When  —  Kit'  it  follows  that  "  Thus  the  duration  until  the  first  tran¬ 

sition  is  U  However,  being  a  constant  \  is  unchanged  hy  the  perturbation,  so 

*By  the  interval  constant  rule,  (he  self  reference  rule  an  !  the  multiplication  rule, 
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Figure  12  Initial  Behavior  Corresponding  to  Path  (QS1.  QS2.  QS3) 
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V'||(0  n;  the  duration  rule  uses  this  fact  to  conclude  that  the  duration  ||.  Hence 
the  conflict  The  perturbation  causes  heat  to  reach  its  transition  quicker,  but 
position  is  unaffected  and  will  transition  at  the  same  time 

Behavioral  inconsistencies  are  located  by  stepping  through  the  transition  in¬ 
tervals  from  earliest  on,  and  checking  the  RC  values  for  the  interval’s  duration 
Section  4.1.1  explains  how  to  find  all  behaviors  that  avoid  this  single  contradic¬ 
tion  while  obeying  the  initial  constraints.  Section  4.1  2  provides  heuristics  for 
eliminating  inferior  paths  Finally,  section  4  1.3  shows  how  to  check  if  the  new 
behavior  is  globally  consistent,  riot  just  a  fix  to  the  first  contradiction.  Note  that 
all  of  these  techniques  depend  <>n  the  DQ  inference  rules  which  are  incomplete 
As  a  result,  while  most  inconsistencies  are  detected,  it  is  not  guaranteed  that 
all  inconsistencies  can  be  found 

4.1.1  Finding  Consistent  Alternatives 

A  simple  observation  about  the  inference  rules  «.f  section  3  forms  the  foundation 
for  the  contradiction  resolution  method  only  the  duration  rule  can  generate  an 
RC  value  for  a  time  duration  Therefore,  the  contradiction  must  be  caused 
by  two  (or  more)  firings  of  the  duration  rule  for  the  same  interval.  What  dis¬ 
tinguishes  these  firings  are  the  different  perspective  parameters  used  in  each 
application  of  the  rule. 

In  the  heat  exchanger  example,  the  two  perspective  parameters  are  Q  and 
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A  In  the  initial  behavior,  they  reached  trans.HK  ns  in  the  same  state,  QS3  But 
to  achieve  consistency  with  t tie  pert urbat i‘  n.  we  must  find  a  behavior  where 
they  reach  transitions  independent iy  this  means  finding  a  path  which  starts 
with  QS1  and  QS2,  arid  passes  through  a  sibling  of  QS3  The  answer,  of  course, 
is  the  path  (QS1  QS2  Q’vV  QSt;>,  QS7i  as  shown  m  figure  13  This  path 
illustrates  the  general  case  A  node  representing  the  qualitative  state  at  a  time 
point  (QS3)  is  replaced  by  three  states:  two  at  time  points  (QS5.  QS7)  and  one 
for  the  interval  connecting  them  iQSbi  tor  the  purpose  of  discussion,  l  shall 
-all  QS3  t tie  FRAGMENTING  point  and  the  two  time-point  states  which  define 
our  objective,  the  PREPOINT  and  POSTPOINT  respectively 9 


Figure  13:  The  Behavior  Corresponding  to  Path  (QS1,  QS2.  QS5,  QS6,  QS7) 


The  problem,  then,  is  to  search  the  state  tree  among  the  siblings  of  the 
fragmenting  point  to  find  the  pre-  and  postpoint  states.  We  know  that  Q  and 
A'  must  reach  transitions  in  different  states  but  which  should  reach  its  transition 
first?  Consider  the  two  duration  RC  values  which  cause  the  contradiction.  Since 
||  specifies  earlier  termination  than  j|  Q,  the  perspective  parameter  for  the  firing 
which  produced  the  fj.  value,  will  reach  its  transition  first  This  means  that  the 


Actually,  this  discussion  assumes  a  simplified  version  of  the  general  problem.  I  assume 
that  the  contradiction  is  caused  by  only  two  firings  of  the  duration  mle,  and  I  assume  that  at 
the  contradiction  can  be  resolved  by  the  addition  of  a  single  new  transition.  The  general  case 
is  a  straightforward  extension.  If  the  QS2  interval  had  three  conflicting  duration  RC  values. 
U,  ||,  and  fl,  then  QS3  could  split  into  five  states  three  for  time  points  and  two  connecting 
intervals.  If  multiple  rule  firings  are  allowed  for  each  RC  value,  then  correspondingly  more 
paths  are  possible. 


preprint  will  have  Q  —  li  an. I  n<>t  .V  =  I'  ^ 1 1 1.- <  QS5  is  tin-  only  stal<-  t . .  nm  t 
this  requirement  so  QS5  is  tin-  |-r-  p.-int  l  or  tin-  In  it  rxiliang-  r  exampl.  .  tins 
state  uniquely  defines  tin  u.-w  I  •«»! ,  ;i\  ior.  b.vau>e  i.ek  ..tie  path  me  I  it*  I*  >  . 

hence  QS7  must  he  tin-  post  point  In  a  inor--  complex  example.  however.  tier- 
could  he  several  candidates  lor  pr-  point  atnl  multiple  behaviors  passing  through 
each  one.  The  following  conditions  further  restrict  the  possibilities. 

•  All  parameters  that  reach  transitions  m  the  prepoint,  must  have  reaclu  l 
transitions  in  the  fragmenting  point 

•  All  parameters  that  reach  transitions  in  the  post  point.  miM  have  reached 
transitions  in  the  fragmenting  point 

•  All  parameters  that  reached  transitions  in  the  fragmenting  point  itnt-t 
reach  transitions  in  either  the  prepoint  or  the  postpoint,  hut  tu>t  both 

While  these  conditions  are  loyal  to  the  implicit  constraints  result  mg  from  the 
initial  selection  of  behavior,  they  are  unfortunately  not  sufficient  to  guarantee 
a  unique  alternate  behavior.  The  next  section  explains  a  heuristic  that  will 
guarantee  a  unique  behavior  hut  not  necessarily  one  that  obeys  all  in/plicit 
constraints. 

4.1.2  A  Heuristic  For  Eliminating  Behaviors 

The  conditions  listed  above  produce  a  unique  behavior  except  in  cases  where 
additional  parameters  besides  Q  and  X  reach  transitions  in  the  fragmenting 
point.  When  extra  parameters  reach  transitions  in  the  fragmenting  point,  one 
must  choose  where  they  should  transition  —  in  the  prepoint  or  the  postpoint. 
The  following  cases  result: 

•  The  parameter  could  be  causally  connected  to  either  Q  or  X.  In  fact, 
this  is  the  case  with  the  heat  exchanger:  F  transitions  to  zero  in  QS.'f, 
How  did  we  know  that  F  should  reach  its  transition  in  the  prepoint  ratlu-r 
than  the  postpoint?  We  didn’t  even  need  to  consider  the  question.  By 
constructing  the  state  tree,  QSIM  already  handled  the  problem  for  use. 
It  recognized  that  F  must  transition  whenever  Q  transitioned:  thus  the 
state  tree  contains  only  this  possibility.  Since  the  topological  consistency 
code  searches  the  state  tree,  it  automatically  benefits  from  QSIM’s  wotk. 

•  There  could  he  additional  RC  information  about  the  parameter.  For  sim¬ 
plicity.  t his  case  was  not  discussed  above,  hut  suppose  that  the  duration 
rule  had  fired  three  times  with  Q ,  A’,  and  S  as  perspectives.  If  .S’  and  X 
both  caused  the  duration  rule  to  deduce  an  RC  value  of  ||.  then  both  .V 
and  S  should  reach  transitions  in  the  postpoint.  Unfortunately,  other  R( ' 
values  complicate  the  analysis.  If  three  different  RC  values  result  from 
the  three  firings,  then  the  fragmenting  point  will  split  into  five  states. 
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The  details  are  messy,  hut  the  concepts  for  resolution  are  similar  to  those 
described  above. 

•  The  parameter  could  be  independent  of  the  perturbation  without  the  in¬ 
ference  rules  deducing  this.  As  with  the  previous  case  (where  the  inde¬ 
pendence,  duration  |j.  was  deduced)  the  parameter  should  transition  in 
whichever  state  has  duration  |j. 

•  The  perturbation  could  change  the  parameter's  transition  time  without 
the  inference  rules  deducing  this.  The  correct  behavior  is  not  predictable 
since  the  change  in  duration  is  not  known 

Since  there  is  no  way  to  correctly  handle  the  last  case,  a  reasonable  heuristic 
is  to  assume  t  hat  it  never  happens.  This  corresponds  to  Occam's  Razor.  Assume 
that  unless  the  duration  rule  says  otherwise,  the  perturbation  does  not  change 
the  transition  time  of  any  parameters.  Thus  if  the  heat  exchanger  example 
had  an  extra  parameter,  ,S\  which  reached  a  transition  in  QS.d,  then  we  should 
assume  that  .S’  transitions  with  A"  in  QS7. 

4.1.3  Ensuring  Global  Consistency 

Using  the  heuristics,  the  algorithm  described  above  is  guaranteed  to  find  a 
unique  postpoint.  But  there  may  be  several  state  tree  paths  that  pass  through 
this  post  point.  To  locate  a  single  new  behavior,  the  program  must  step  through 
the  original  behavior  from  the  fragmenting  point  onwards.  Every  time  a  branch 
in  the  tree  is  taken,  the  corresponding  descendant  of  the  postpoint  should  be 
selected  as  well.  When  the  original  behavior  reaches  a  leaf,  a  unique  new  be¬ 
havior  will  result.  Unfortunately,  there  are  two  reasons  why  processing  must 
continue. 

•  Many  RC  values  must  be  recomputed.  Because  the  RO  values  refer  to 
transition  points  and  intervals,  all  values  from  the  fragmenting  interval 
onward  will  be  incorrect,  rhis  isn't  very  surprising;  after  all,  we  started 
with  conflicting  duration  RC  values  in  the  first  place  Given  the  new 
behavior,  the  inference  rubs  of  section  3  must  be  rerun  to  generate  a 
consistent  set  of  R(  ’  values. 

•  What  if  these  rules  generate  a  new  contradiction?  There  is  no  guarantee 
that  the  new  behavior  is  topologically  sound  However,  if  conflicting  dura¬ 
tion  RC  values  err  gem-rated  for  an  interval,  that  interval  must  occur  after 
any  interval  which  caused  a  previous  conflict  T  hus  each  cycle  of  inference 
rules  and  topology  resolution  guarantees  that  the  time  of  first  inconsis¬ 
tency  increases.  Since  all  behaviors  are  finite,  the  cycle  must  eventually 
terminate. 
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•  SPLITTING 

Splitting  happens  when  the  duration  rule1"  deduces  a  single  RC  value 
of  ft  or  1J.  fro. n  two  different  perspectives  For  example,  consider  the  heat 
exchange  with  an  initial  behavior  of  (QS 1 ,  QS2.  QS.'l)  and  the  perturbation 
of  —  A'ft  and  l 'ff .  Since  thermal  conductivity  is  higher,  equilibrium  will 
occur  sooner,  but  since  the  oil  is  moving  faster,  it  will  get  out  quicker.  Both 
parameters  lead  the  duration  rule  to  conclude  QS3  will  occur  quicker;  thus 
the  initial  behavior  is  consistent  But  so  is  every  behavior.  If  the  thermal 
conductivity  is  much  higher  and  velocity  is  only  a  little  higher,  then  the 
behavior  (QS1.  QS2.  QS.').  QS(>,  QS7)  will  result  If  l'  was  increased 
more  than  —  K  then  the  path  (QSl  QK2.  QS  1)  would  result.  Because 
the  perturbation  was  specified  in  qualitative  terms,  there  isn’t  enough 
information  to  resolve  the  ambiguity  and  CA  must  return  all  possible 
splits  of  the  two  parameters  Q  and  X 

Like  the  techniques  of  section  1.1.  my  methods  for  finding  other  consistent 
behaviors  are  dependent  on  the  DQ  inference  rules.  As  a  result  they  are  neither 
complete  nor  sound.  For  example,  suppose  the  duration  of  an  interval  was  (f. 
but  the  duration  rule  had  not  deduced  this  fact.  Then  compaction  would  not 
be  considered  and  a  possibly  consistent  behavior  would  not  be  considered  Sim¬ 
ilarly,  one  of  the  techniques  could  suggest  a  behavior  which  appears  consistent 
only  because  the  I)Q  rules  were  inadequate  to  expose  a  contradiction. 


°S|)litting  is  the  only  ra.sc  (bat  anal\7cs  just ifirat i..ns  ami  depends  on  the  fart  that  the 
duration  rule  is  the  only  way  to  generate  a  duration  RC  value.  Comparting,  stalling,  and 
kick-starting  only  require  the  RC  value  and  access  to  the  state  tree. 


5  Related  Work 


Although  comparative  analysis  questions  have  long  been  important  topics  in  t  In¬ 
fields  of  engineering  and  mathematics,  little  work  has  been  done  on  comparative 
analysis  in  the  artificial  intelligence  community 

5.1  Sensitivity  Analysis 

Sensitivity  analysis  is  a  common  engineering  technique  for  calculating  the  effect 
on  system  performance  due  to  variations  in  system  parameters.  In  other  words, 
comparative  analysis  is  a  qualitative  version  of  sensitivity  analysis.  The  sensi¬ 
tivity  of  a  quantity  T,  to  perturbations  in  a  parameter  A',  is  defined  [1]  as  tin- 
prod  net : 

OT  X 
ct.X  T 

Because  of  its  important  application  to  design,  considerable  work  has  been  done 
on  efficient  methods  for  calculating  sensitivities.  Approaches  include  numeri¬ 
cal  and  symbolic  differentiation,  construction  of  an  incremental  network,  and 
analysis  of  an  adjoint  network  [1], 

Compared  to  DQ  analysis,  these  methods  have  a  major  advantage — they 
generate  a  quantitative  value  for  sensitivity.  But  sensitivity  analysis  has  two 
limitations:  it  does  not  generate  explanations,  and  it  requires  an  explicit  equa¬ 
tion  for  tiie  desired  quantity  T.  Thus  sensitivity  analysis  could  not  solve  the 
spring/block  problem  until  the  human  modeler  provided  a  formula  for  period. 

The  technique  of  comparative  statics  [17,12],  long  used  in  economics  to  com¬ 
pare  two  different  equilibrium  behaviors,  suffers  from  the  same  limitation  It 
requires  explicit  formulas  for  the  partial  derivatives  in  question 


5.2  Partial  Derivatives 


Since  the  RC  notation  expresses  how  a  parameter  changes  given  an  initial  per¬ 
turbation,  it  is  natural  to  ask  about  its  relationship  to  the  standard  mathemat¬ 
ical  tools  for  expressing  relative  change:  partial  derivatives.  In  the  following 
proposition  it  is  handy  to  think  of  parameter  C  as  the  cause,  and  E  as  an  effect. 

Proposition  15  If  C ff0  and  all  other  independent  and  boundary  condition  pa¬ 
rameters  have  an  RC  value  of  ||0  and  E- lj0  then 


dE  n 

dC  <  ° 


at  time  zero. 


This  statement  can  be  extended  to  any  transition,  y,,  by  normalizing  with 
respect  to  time.  While  the  relationship  between  RC  values  and  partial  derivative 
is  straightforward  for  values  at  transition  points,  the  connection  is  more  subtle 
for  interval  RC  values  because  of  the  presence  of  perspect  ives. 


5.3  QP  Theory 

In  his  treatise  oil  Qualitative  Process  theory  1  orbus  discussed  difrereiit  ini  quali¬ 
tative  analysis  [7,  pages  159-lfil],  I>ut  attempted  no  implementation  He  delined 
quantities  g!  greater  than  < over  an  interval,  i.  if  for  all  instants  in  the  interval. 
q\  >  q 2  measured  at  that  instant  C lifort unately.  this  definition  has  several 
problems.  Since  the  quantification  is  over  a  single  interval  of  time,  it  is  impos¬ 
sible  to  make  comparisons  of  systems  whose  tune  behavior  changes  as  a  result 
of  a  perturbation.  Thus  his  attempt  to  formalize  "distance  equals  rate  times 
duration”  in  predicate  calculus  is  severely  limited  Rates  can  only  be  compared 
if  the  duration  of  an  interval  is  unchanged' 

But  even  if  the  quantification  was  correct,  time-wise  comparison  is  almost 
never  a  useful  one  to  make.  In  the  spring/ block  case,  for  example,  it  simply  isn't 
the  case  that  the  heavy-block  is  always  moving  slower  than  the  small-block;  the 
periods  get  out  of  phase.  The  key  to  solving  these  problems  is  in  the  use  of 
perspectives,  discussed  in  this  paper.  The  comparison  on  velocity  (necessary  to 
predict  that  the  period  lengthens)  is  valid  only  from  the  perspective  of  position 

5.4  Temporal  Representation 

QSIM  [15]  is  an  efficient,  easy  to  use  simulator  that  has  significantly  sped  the 
development  of  both  my  comparative  analysis  theory  and  the  CA  implementa¬ 
tion.  However,  QSIM  has  defects;  its  weak  temporal  representation  is  a  major 
problem. 

As  explained  by  Hayes  in  [11],  systems  which  represent  behaviors  as  a  se¬ 
quence  of  states  force  a  total  ordering  on  events.  Because  qualitative  reasoning 
is  often  unable  to  unambiguously  determine  an  order,  the  behavior  must  branch 
to  consider  multiple  possibilities  If  events  interact,  then  the  various  branches 
often  have  interesting  qualitative  differences  But  frequently,  the  alternate  be¬ 
haviors  are  equivalent  and  just  complicate  reasoning  and  consume  processing 
resources 

To  combat  this  problem.  Williams  introduced  the  notion  of  concise  episodes 
[23],  and  has  devised  an  efficient  simulator  (railed  a  Temporal  Constraint  Prop¬ 
agator)  to  manipulate  them  J list  as  qualitai  ive  simulators  using  Williams'  tem¬ 
poral  representation  would  improve  on  QSIM.  comparative  analysis  programs 
would  have  several  advantages  over  CA  Williams  is  building  such  a  system  for 
use  in  automated  design  [25] 

•  The  propositions  of  section  .'(  would  still  be  true,  and  could  be  encoded 
more  easily.  CA  requires  explicit  rules  for  composing  durations  over  inter¬ 
vals  (e  g.,  if  t)URATIONff(Q  ,  ,  and  I)l  RATION/)  ,  J,  then  nt'RATIONf[(0  2|). 
These  computationally  expensive  rules  would  be  subsumed  by  the  tempo¬ 
ral  constraint  propagator 


The  search  for  topologically  distinct  behaviors  (section  4)  would  be  simpli¬ 
fied  because  the  space  would  be  smaller.  By  eliminating  irrelevant  order 
distinctions,  the  number  of  different  behaviors  would  be  smaller  Only 
if  it  was  qualitatively  interesting  would  there  be  any  need  to  consider  a 
behavior  in  which  two  parameters  reach  transitions  simultaneously 


6  Future  Directions 


This  paper  discussed  the  problem  ■>!' romparat ive  analysis.  the  task  oi  explaining 
how  a  system  will  react  to  perturbations,  and  why.  Multiple  perspectives,  which 
can  be  used  to  reparaineterize  system  parameters,  lead  to  a  powerful  definition  of 
relative  change.  DQ  analysis  solves  comparative  analysis  problems  by  applying 
inference  rules  to  the  initial  perturbation  of  a  system  A  trace  of  the  rules  used 
in  solving  a  problem  can  be  easily  translated  into  an  intuitive  explanation  of 
the  answer  Since  the  rules  have  been  proven  sound.  DQ  analysis  is  guaranteed 
to  produce  only  correct  explanations.  A  computer  programs,  CA.  implements 
the  theory  of  DQ  analysis  and  correctly  solves  over  twenty  comparative  analysis 
problems  including  those  that  change  the  order  of  transitions  in  the  behavior. 

Despite  the  success  of  DQ  analysis,  several  areas  for  future  research  beckon. 

•  DQ  analysis  is  incomplete.  Although  CA  is  guaranteed  to  terminate,  it 
doesn't  always  deduce  an  RC  value  Fortunately,  there  are  other  tech¬ 
niques  for  solving  comparative  analysis  problems  Exaggeration,  for  ex¬ 
ample,  saves  many  problems  with  a  completely  different  style  of  reasoning 
[20],  Although  it  is  believed  that  exaggeration  is  also  incomplete,  initial 
results  suggest  that  exaggeration  can  solve  several  problems  which  DQ 
analysis  cannot  [22] 

•  Certain  comparative  analysis  questions  have  no  answer  For  example, 
“What  happens  to  the  period  of  oscillation  of  a  spring/block  if  both  the 
mass  and  spring  constant  are  increased?”  It  would  be  nice  if  C‘A  could 
recognize  that  there  was  no  answer  to  this  question  instead  of  simply- 
saying  that  it  can't  find  an  answer. 

•  Any  analytic  technique  is  only  as  good  as  the  model  on  which  it  works. 
Currently,  humans  construct  models  and  computers  are  only  used  in  ana¬ 
lyzing  them.  This  imparts  fragility  to  the  process  For  example,  consider 
the  structural  description  of  the  oscillating  spring/block  example  (section 
3.8).  Suppose  that  the  initial  situation  had  X  =  0,  V  set  to  some  max¬ 
imum  value  and  the  perturbation  was  Cff  Although  the  DQ  inference 
rules  can  deduce  facts  like  “the  maximum  displacement  will  increase,”  the 
topological  analyzer  is  unable  to  recognize  the  possibility  that  the  spring 
will  break. 

The  cause  is  a  simplistic  model  Hooke's  law  precludes  the  possibility  of 
a  broken  spring.  To  achieve  greater  robustness  in  qualitative  analysis  in 
general  and  comparative  analysis  in  specific,  modeling  must  be  treated  ex¬ 
plicitly.  By  incorporating  ontological  assumptions  into  process  definitions, 
QP  theory  [7]  has  made  progress  here,  but  further  research  is  necessary  to 
address  the  questions  of  reasoning  with  multiple  models,  dynamic  model 
creation,  and  the  evaluation  of  a  model  in  t  he  context  of  a  specific  problem 
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A  A  Useful  Example 

This  section  constructs  an  example  which  serves  both  as  a  counter-example  for  a 
generalized  version  of  the  derivative  rule  and  as  the  proof  of  the  non-exhaust  ive 
proposition.  Suppose  that  .1,  V.  and  A  are  parameters  such  that  .1  =  \ 

V  =  jj-  A',  and  A’  is  a  covering  perspective  over  -;.,+  ] ),  The  derivative  rule 
(proposition  3)  showed  that  if  .4  and  V'  are  positive  over  the  interval  -;,  +  1  j 
and  if  -d'-fr,,  and  A^ll+ ,,,  then  1+1). 

Unfortunately,  the  derivative  rule  is  not  true  for  arbitrary  perspectives.  The 
following  abberation  should  convince  you  of  this.  I  show  thr<'e  parameters.  U.  .1 
and  P  such  that  .4  =  4-  V’  and  P  is  a  covering  perspective  over  (*•„.  "h  )•  Yet  al¬ 
though  ,4||^  j  j  the  parameter  1'  has  no  consistent  behavior  from  the  perspective 

of  P  During  part  of  the  interval  V and  during  part  l 

Here  are  the  details.  Over  the  absolute  time  interval  (0.  1)  define: 

no  =  kr- 

A(t)  =  &H0  =  <  -4  ” 1  (a )  =  n 
P(t)=A(t)  =  t  P~l(p)  =  P 

v(t)  =  k* 3 

A(t)=  &V(t)=  t2  .4-‘(a)  =  ui 
P(t)  =  A(t)  =  l2  P~'(p)  =  pl 


Note  that  /’(0)  =  P(0)  =  0  and  /’( 1 )  =  /*(  1 )  =  1  and  P  is  strictly  monotonic, 
so  P  is  a  valid  perspective  over  this  interval.  Since  P  —  A  and  P  =  .4  the  self 
reference  theorem  shows  that  .4||^  lr  So  what  does  V  do  from  the  perspective 
of  P?  Consider  p  =  C 

t  i,!  =  ^.  a,, i  =  =  l 


Now  let  p 


v(r-'(p))=  -p- 


j|y-  and  V(P  '(/>))  = 


So  for  a  small  value  of  p  the  corresponding  f  is  larger  than  r  but  for  larger  p 
t he  situation  is  reversed.  Thus  it  is  neither  t he  case  that  Uff^  , ,  nor  l’||  „  n  nor 

PV<'n  •hough  AW^ 
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